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ABSTRACT 


All Navy aircraft are required to accommodate, 
anthropometrically, ninety percent of the user population. 
Some designs have been criticized for their low accommo- 
dations but those accommodations have never been quantified. 
The purpose of this thesis was to quantify the accommoda- 
tion, by each type of operational Naval aircraft, of 
populations of Naval aviation personnel of 1964, 1969, and 
1975. The Computerized Accommodation Percentage Evaluation 
(CAPE) model was used to generate data points since only 
summary statistics were available for two of the popula- 
tions. Each subject of every population was checked against 
the requirements of the design specification, and against 
the limitations of each aircraft. All aircraft were found 
to accommodate more than ninety percent of the 1975 popu- 
lation. Time related changes in the populations were noted 
and unexplained inconsistencies in the data were discovered. 
Possible sources of error were discussed and potential 


solutions proposed. 
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iP CRObDUCT TON 


The Human Engineering Design Criteria for Military 
Systems, Equipment and Facilities [MIL-STD-1472B] states 
the requirement for anthropometric accommodation of 
equipment users: 

Design and sizing shall insure accommodation, com- 
patibility, and maintainability by at least 90 
percent of the user population. Generally, design 
limits shall be based upon a range from the 5th 
percentile to the 95th percentile values for 
critical body dimensions. 

For the special case of Naval Aircraft cockpits, the 
accommodation criterion was enlarged in 1973 to include the 
central 95 percent of the user population [CNAVAIR SD-24K, 
1973]. The anthropometric description of the population is 
based on data collected in 1964. Although a future survey, 
scheduled for 1981, will provide a more current description 
of the anthropometric features of Naval aviation personnel, 
aircraft now being designed for production and use well into 
the future are being designed to accommodate the 1964 popu- 
lation. The F-18, projected to become operational in 1982 
will be built to accommodate the pilots of 1964. Other 
designs, if they have advanced much beyond the concept 


stage, are committed to the same aircrew accommodation 


standards. 


ipl 





The anthropometric characteristics of the population 
are changing but the design specification is not. That the 
population is changing is reflected in efforts by the 
National Aeronautics and Space Administration and by the 
United States Air Force to predict the typical pilot size 
of the future. A 1978 study conducted by NASA concluded 
that stature (height) increases among pilot and potential 
astronaut males could be expected to amount to about 8 mm 
(1/3 inch) per decade [NASA 1024 Vol. I, LOW Su LG 
U.S. Air Force predicted an increase in height of 0.418 
inch per decade Roebuck, Kroemer, and Thompson, 1975 
These changes, since 1964 and into the future, affect 
pilot accommodation in operational Naval aircraft, but the 
exact nature of the effect is not known. 

In 1977 the Navy implemented the Anthropometric 
Compatibility Assignment Program [OPNAVINST 3710.36], a 
system aimed at preventing the assignment of Naval Aviators 
and Naval Flight Officers to aircraft which are incompatible 
with their anthropometric features. Under the provisions of 
this program, pilots' and NFOs' anthropometric measurements 
are taken when they enter flight training ({BUMEDINST 3710.1, 
1977]. Values for sitting height, functional reach, 
buttock-knee length, and leg length (for descriptions of 
anthropometric terms, see Appendix A) are encoded and 


become part of their permanent records. For example, a 
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man coded 

G74 2°33 
would have a sitting height measurement between 38.95 and 
Bens inches, a functional reach of 30.5 to 30.9 inches, 
buttock-kKnee length of 23.0 to 23.9 inches and a leg length 
oeeeO.0 to 42.9 inches. 

The cockpits of each operational aircraft type are 
measured to determine the limits of pilot size compatible 
with the cockpit ENAVAIRINST 3710.9, 1979]. The maximum 
and minimum acceptable dimensions for sitting height, 
functional reach, buttock-knee length, and leg length are 
determined for each cockpit. These maxima and minima are 
then converted to four coded descriptors, which signify 
dimensions of exclusion for pilots. For example, an A-4E 
ds coded: 

089 9 789 0 
AS with the personal anthropometric codes, the first of the 
four groups of numbers corresponds to sitting height. The 
second group corresponds to functional reach; the third to 
buttock-kKnee length, and the fourth to leg length. The 
code 089 9 789 O indicates that a pilot with a sitting 
height code O, 8, or 9 would not be eligible for assignment 
to an A-4E. Similarly any pilot with a functional reach 
code 9, buttock-knee length code 7, 8, or 9, ora leg 
length code O would be excluded. 


The Chief of Naval Aviation Training (CNATRA) is 
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required to utilize these codes early in each aviation 
student's training to exclude the student from training 
sequences leading to aircraft with which he is not com- 
patible [CNATRAINST 13520.1, 1980]. The Navy Military 
Personnel Command uses the anthropometric codes in the 
process of assigning pilots and NFOs to aircraft squadrons. 

Currently, a similar coding scheme is under consider- 
ation for use internationally among free world nations 
(Clauser, 1980]. The need for such a program was indi- 
cated by difficulties encountered by some foreign countries 
operating aircraft acquired from the United States and 
other foreign sources. 

In Figure 1, it is apparent that the height of U.S. 
Air Force flight personnel is generally greater than that 
of other countries. The fifth percentile in U.S. height 
is greater than the SOth percentile of the samples from 
Japan, Thailand, and Vietnam. Figure 2 shows a similar 


melationship for sitting height. 
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Peanoui ea. 
Heights of Seven Military Populations 
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Figure 2 


Sitting Heights of Seven Military Populations 
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In this case the Sth percentile is larger than the 50th 
percentile of Thailand and Vietnam. Aircraft designed to 
accommodate people as small as the Sth percentile of the U.S. 


population can be expected to exclude about half of certain 
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foreign populations and as much as 75 percent of the 
Vietnam aviation personnel [Kennedy, 1973]. 

Although design criteria require accommodation of the 
central ninety percent of the Naval aviation population, 
the existing population is imperfectly described by the 1964 
sample. Changes in the selection process, trends in the 
United States general population were cited as factors con- 
tributing to the differences between the 1964 and the 1969 
populations [Moroney, Kennedy, Gifford, and Provost, 1971]. 
In addition, designing to accommodate the fifth and ninety- 
fifth percentiles in each variable does not yield an 
accommodation of ninety percent of the total population. 
Since exclusion results from any one measurement being out- 
Side limits, accommodation is as much a function of the 
correlation between measurements as it is a function of 
each measurement. The correlation between measurements is 
appreciably less than 1.0 and accommodation is less than 
that indicated by any univariate estimate. When twelve 
measurements were considered, the exclusion of subjects who 
were outside the Sth and 95th percentile range in any 
measurement yielded a cumulative exclusion of 48 percent 
[Moroney and Smith, 1972). As the existing population is 
altered further by expanding the acceptable range to 
include a greater percentage of the female population, user 
accommodation will be adversely affected. For example, a 


generalized cockpit built to male anthropometric standards 
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excluded 22 percent of the male user population. The same 
cockpit excluded 90 percent of the females in the U.S. Air 
Force [Ketcham-Weidl and Bittner, 1976]. 

Before considering future changes in the population, 
the current anthropometric exclusion situation should be 
determined. However, the proportion of the Naval aviation 
population excluded by the existing design specification 
has not been quantitatively determined. Similarly the 
proportion excluded by each specific aircraft design is 
unknown. Some aircraft can be expected to exclude more 
than the ten percent designed out by the specification. 
Other aircraft will exclude less than the maximum allowable 
ten percent by virtue of design allowance or perhaps by 
chance. 

The methods employed in measuring aircrew personnel 
have been found to be rather unreliable [Moroney, et. al., 
1971]. These measurements, however, form the basis of the 
Anthropometric Compatibility Assignment Program and their 
accuracy is critical since they constitute a determining 
maectOr in aneofficer's career path starting point. Addition- 
ally, these same data will be used to help identify (or 
eliminate) subsequent aviation duty assignments. Consider- 
able attention has been placed on determining the proper 
use of the data but the actual collection of the data 
appears to be an area where reliability could be improved. 


Currently, there is no systematic method of checking 
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anthropometric measurements. Each student is measured once 
without any attempt at remeasurement or checking of the 
accuracy of the measurements. If a measurement is recorded 
in error, the error becomes part of the student's record, 
undetected until some independent event brings the error to 
light. Errors in measurement of as much as three inches 
have been found [O'Leary, 1980] subsequent to the erroneous 
measurement being used to determine anthropometric suit- 
ability and assignment of a student to an aircraft type. 
Some students have completed flight training before anthro- 
pometric measurement errors were detected, resulting in 
reassignment and expensive retraining. 

The purpose of this thesis is to quantify the amount 
of incompatibility between aircraft and aircrewmen and to 
propose some alternatives to improve the reliability of the 
personal anthropometric measurements now being taken. By 
examining samples drawn from three Naval aviation popula- 
tions, each of a different year, proportions of those 
populations excluded from Naval aircraft will be determined. 
The exclusions resulting from changes in the population will 
be quantified and the specific proportions excluded from 
each aircraft type, as a function of sample source, will be 
calculated. Those results will lead to a discussion of some 
inconsistent patterns in the data and recommended changes in 


the data collection procedure now used by the Navy. 


Ve, 





Chapter 2, DATA SOURCES, describes the three data sets 
analyzed. For each of the three samples: the subjects were 
described, previous statistical analyses discussed, and 
methods of data collection compared. Chapter 3, DATA 
PREPARATION, describes the preparation of each data set 
to facilitate analysis. Preparations consisted primarily 
of the use of a segment of the Computerized Accommodation 
Percentage Evaluation (CAPE) model to generate data points 
from a statistical description of the population. Chapter 
4, ANALYSIS OF DATA, describes the screening of each data 
set to determine the proportion of the sample excluded. 
Secreening involves the elimination of any subject with an 
anthropometric measurement outside prescribed limits. The 
proportions excluded from each sample are then compared. 
Observations on the apparent shortcomings in data collection 
methods, errors in measurements, and unexplained inconsis- 
tencies were included. Chapter 5, SUMMARY AND CONCLUSIONS, 
includes inferences concerning interpretations of anthropo- 
metric data and recommendations in the area of data 


collection methods for use by the Navy. 
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Lie Ad ee SOURCKHS 


The data to be analyzed consisted of samples from three 
populations. Samples of Naval aviation personnel were taken 
in 1964, 1969, and 1975. Differences between populations 
other than the year of the sample exist and they will be 


discussed below. 


A. THE 1964 SURVEY DATA 

In 1964, a survey was conducted to determine the anthro- 
pometric chhamagiseren areas of personnel flying in Naval 
aircraft. The survey was conducted in a manner to yield 
results sufficiently precise to be utilized subsequently for 
design of aircraft. The Bureau of Medicine and Surgery, 
tasked to collect data for use by the Bureau of Weapons, 
sent a trained anthropometric measuring team to ten different 
Naval and Marine Corps air stations to take measurements. 

Included in the sample were pilots, Naval Aviation 
Observers, bombardier-navigators, radar observers, flight 
surgeons who were designated Naval Aviators, and enlisted 
personnel who were permanently designated aeronautical 
personnel and who were still eligible for assignment to 
subsequent duty involving flying. Naval Aviators included 
both U.S. Navy and the U.S. Marine Corps flyers. 


The subjects in this survey comprised approximately ten 
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percent of the Naval aviation population on active duty in 
1964. The sample included personnel of various ages, ranks, 
ethnic backgrounds and sources (Naval Academy, NROTC, etc. ). 
In that respect the sample undoubtedly provided an accurate 
representation of the Naval aviation population. The com- 
position of the sample is somewhat different from the later 
samples in the respect that enlisted men and older men were 
surveyed in 1964 but not in the later samples. The relia- 
bility of the measurements received special attention and 
was enhanced by the repeat measurement of randomly selected 
subjects, insuring rapid detection of measurement errors and 
encouraging the data collectors to be careful in their task. 
Each measurement was taken by personnel specifically trained 
for the job, using a technique of measurement recognized as 
the most accurate available. 

Published results of the survey include means and 
standard deviations, as well as percentile listings for each 
variable [Gifford, Provost, and Lazo, 1965]. Subsequent 
analysis of the data led to publication of correlations 
between pairs of variables (Moroney, et. al., 1971]. 

The results of this survey were adopted by the Navy as 
the basis for aircraft design specifications. MIL-STD-1333A, 
Military Standard Aircrew Station Geometry for Military 
Aircraft, includes NAEC-ACEL-553 as part of the standard, 
thus specifying the 1964 data as the guideline for aircraft 


being built today (NAEC-ACEL-553 reported the results of the 
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1964 survey). In that sense, these results are recognized 
as the standard and they have become the yardstick against 
which other data can be compared. 

This data set has two shortcomings with respect to this 
thesis, however. No measurements of leg length (or buttock- 
leg length) were taken. The components of leg length 
(buttock-knee length, knee height sitting and popliteal 
height sitting) were recorded but leg length had not yet 
been recognized as a critical selection dimension. Also, 
the results of data analysis are available, but the raw 
data are not. The importance of these two limitations will 


become apparent when procedures are discussed. 


By THE 1969 DATA 

Subsequent to the introduction of an integrated anthro- 
pometric measuring device in 1964, the Navy required that 
the device be used to measure all aviation training candi- 
dates during their aviation training entrance physical 
examination, which was administered at the School of Aviation 
Medicine (now the Naval Aerospace Medical Institute). The 
collection of anthropometric data was thus institutionalized 
and became a routine process. Eight measurements were taken 
on each entering Student Naval Aviator and Student Naval 


Flight Officer. These measurements were: 


weight trunk height 

height buttock-Knee length 
sitting height buttock-leg length 
Shoulder width functional reach 
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Between January 1966 and August 1969, 6,534 students 
were measured. Their measurements comprise the data set 
herein called the 1969 data. The description of this data 
set was published in 1971 CLMoroney, et. al., 1971] and 
includes the mean and standard deviation of six measurements 
(shoulder height and functional reach were deleted from the 
analysis), percentile listings for those measurements and 
correlations between variables. 

No information on functional reach or shoulder height 
was published because of the discovery of measurement tech- 
nique errors resulting in measurements uncorrectible by 
any means short of repeat measurement. This deviation from 
ideal procedure should not have been altogether surprising 
Since measurements were taken by personnel whose level of 
training did not match that of the 1964 data collectors. 

In addition, measurements were taken and recorded with no 
subsequent checking, leaving procedural errors undetected 
and biasing the measurements of many subjects. 

The size of the 1969 data set is quite large, consisting 
of virtually all of the student aviation input for over three 
years. AS with the 1964 data set, no raw data is available, 
only analytical results. The absence of functional reach 
from the data analysis will present an obstacle to be 


discussed in Chapter 3. 
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See ide 1975 DATA 

Following the publication of the results of the 1969 
data analysis, an improved integrated anthropometric 
measuring device was installed at NAMI and it was subse- 
quently used to measure incoming students. That device is 
still in use today. It is operated by NAMI staff Hospital 
Corpsmen and it was used to measure the same eight dimen- 
Sions as was its predecessor. A sample of 968 students 
who entered training in 1975 comprise the subjects of the 
1975 data. The raw data is available and from an examina- 
tion of the data points thirty-two subjects were found to 
be deficient in some respect and did not meet the entry 
requirements delineated by the Bureau of Medicine and Surgery 
[LMANMED, 1978]. Most of these ineligibilities were based on 
excessive weight. However, all subjects were screened for 
minimum and maximum values of weight, height, sitting height, 
and buttock-leg length, as well as for height-weight rela- 
tionship. The remaining 936 subjects comprise the acceptable 
sample. 

The 1975 data were collected in a manner similar to that 
used for the 1969 data set: a routine execution of estab- 
lished procedures by Navy enlisted personnel. The 1975 
data also manifest data collection deviations similar to 
those of the 1969 data. The number of subjects is smaller 
than either of the other samples but still large enough to 


lend itself to statistical analysis. 
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Till. DATA PREPARATION 


The goal of the analysis was to determine the proportions 
of the 1964, 1969, and 1975 samples excluded by the four 
anthropometric limitations (sitting height, functional 
reach, buttock-knee length, and leg length) required by 
OPNAVINST 3710.36A. This goal was achieved by screening 
each subject four times, once for each of his measurements, 
and eliminating every subject whose measurements were not 
weeeeewithin limits. 

This procedure required Knowledge of the four individual 
measurements for each subject. In the 1975 data set, the 
four measurements (and more) were readily available. For 
the earlier data sets, only means, standard deviations, and 
correlations were available and no individual measurements 
were known. It was necessary, therefor, to generate data 
to simulate those measurements. 

The Computerized Accommodated Percentage Evaluation 
(CAPE) model [—Bittner, 19751] was used to generate the 
required data points. This model required as entering 
arguments the mean and standard deviation for each of the 
four variables, as well as a matrix of correlations between 
variables. The model used a Monte Carlo process to generate 
data points. A program listing is included in Appendix C 


(Appendix D presents the same program modified for use with 
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the IBM 360-67 computer and other machines using Fortran IV). 
The CAPE model has been incorporated in the Computerized 
Accommodated Reach Model and is available at the Naval 
Aviation Development Center, Warminster, Pennsylvania. 

Unfortunately, some of the values required for generating 
the data points were missing from the survey reports. 
Specifically, the 1964 data included no information on leg 
length and the 1969 data omitted functional reach. It was 
necessary to find values in similar populations and substitute 
them for those missing. Taking values from other populations 
required the assumption that concurrent populations, similar 
in cultural, racial, and demographic composition, and having 
passed similar anthropometric entrance requirements, will 
have comparable anthropometric features. That assumption 
was made and surrogate values for the missing means and stan- 
dard deviations were taken from populations which matched 
the descriptions of the populations under analysis (1964 and 
1975), both subjectively and numerically. Inter-variable 
correlations were determined on much the same basis: substi- 
tuting correlations from similar populations. Ketcham-Weidl 
and Bittner [£1976] had substituted correlations from a data 
set of male subjects into a set of female data and exper- 
ienced little loss of precision. 

For the 1964 data, the mean and standard deviation for 
leg length were assumed to be 43.82 and 2.01 inches, respec- 


tively. Those values were established in a 1968 survey of 
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U.S. Army helicopter students ESchane, Littell and Moultree, 
1969]. Of the information available, the population of that 
Survey, better than any other, approximated the 1964 Naval 
aviation population with respect to the factors which Moroney 
[1971] reported as influencing differences between the 1964 
and 1969 Navy data sets. The Sitting height-Leg length cor- 
relation maintained the same value (0.45) for the 1969 and 
1975 data sets so the same value was used for 1964. Func- 
tional reach-Leg length and Buttock knee-Leg length 
correlations (0.639 and 0.860) were extracted from a 1970 
survey of U.S. Army aviators EChurchill, McConville, Laubach, 
and White, 1971], again because of the similarities in 
descriptions of the 1970 Army and 1964 Navy populations. 

For the 1969 data, Functional reach values for USAF 
flight personnel reported by Churchill C1971] were used for 
the mean and standard deviation (31.24 and 1.62 inches). 

The sitting height-functional reach and the buttock knee- 
functional reach correlations were assumed to be the same as 
those found in the 1975 data, 0.44 and 0.63, respectively. 
These values were greater than those of the 1964 data and 
their use, rather than the 1964 values, would influence 
error in the conservative direction, that is, in the direc- 
tion of lesser difference between univariate and multi- 
variate exclusion. The leg length-functional reach corre- 
lation was taken from Churchill's 1970 survey and is the 


same value assumed for the 1964 data. 
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CAPE INPUTS AND OUTPUTS 


Table= i 


FOR THE 1964 DATA 


Means, Standard Deviations and Ranges 


Variable 


poe oins Height 
mean 
std.dev. 
range 


Functional Reach 
mean 
std.dev. 
range 


Buttock-Knee Length 
mean 
std.dev. 
range 


Leg Length? 
mean 
std.dev. 
range 


4Schane, W.P., Littell, 


Observed 
> Dania 


36.28 
ei 
32.19-41.62 


31.51 
1.42 
27.26-36.31 


24.09 
oo 
20.73-27.81 


43.82 
oon 


D.E., and Moultree, 


Generated 
~ Data 


36. 
. 26 
oi 


31. 
.40 
26. 


24 . 
94 
Ze 


43. 
de. 
37. 


C. 


20 


81-40.25 


47 


35-35.38 


10 
0O6-27.11 
82 


99 
71 = oom 


Ge 


Selected Anthropometric Measurements of 1,640 U. Se Army 
Warrant Officer Candidate Flight Trainees, 


Aeromedical Research Laboratory Report 69-2. 
in: Anthropometric Source Book Volume IT: 


Cio A riny. 
as reported 
A Handbook of 


Anthropometric Data, NASA Reference Publication 1024, 


1978, p. 156. 
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Table II 


CAPE INPUTS AND OUTPUTS 
FOR THE 1964 DATA 


Correlations 
Sitting Functional Buttock- Leg 
Height Reach Knee Length Length 


Observed Data 


Sitting Height iO .376 Moo .454 
Functional Reach 1.0 .586 .639) 
Buttock-Knee Length en . 865 
Leg Length 0 


Generated Data 


Sitting Height Pe .410 .374 .446 
Functional Reach a) ooh oo 
Buttock-Knee Length tO oot 
Leg Length a0 


2Same as observed value for 1969 and 1975 samples. 


Dchurchill, E., McConville, J.T., Laubach, L., and 
White, R.M., Anthropometry of U.S. Army Aviators-1970, U.S. 
Army Natick Laboratories Technical Report 72~52-CE, December 
meee, p. 278. 
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Table III 


CAPE INPUTS AND OUTPUTS 
FOR THE 1969 DATA 


Means, Standard Deviations, and Ranges 


Observed Generated 

Variable Data Data 
Sitting Height 

mean 36.76 36.77 

std.dev. ile a eee 

range 31.70-41.0 32.43-40.61 
Functional Reach 

mean 31.24 Si. 20 

std.dev. 162 1.60 

range eee ZU O2=30mrS 
Buttock-Knee Length 

mean 24.54 24.55 

std.dev. 1.26 ie 3 

range 20.40-29.90 20.67-28.31 
Leg Length 

mean 43.86 43.86 

std.dev. 2.08 ZOD 

range SZ 0=00. 10 38.14-50.02 


penmrecnall, £., McConville, J.T., Laubach, L., and 
White, R.M., Anthropometry of U.S. Army Aviators-1970, U.S. 
Army Natick Laboratories Technical Report 72-52-CE, December 
meee, p. G91. 
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Sitting Height 
Functional Reach 
Buttock-Knee Length 
Leg Length 


Sitting Height 
Functional Reach 
Buttock-Knee Length 
Leg Length 


Take a 


CAPE INPUTS AND OUTPUTS 
FOR THE 1969 DATA 


Correlations 
Sole teen Functional Buttock- Leg 
Height Reach Knee-Length Length 
Observed Data 
0 .444 .40 245 
Tec .634 .63g9b 
10) .86 
ee 
Generated Data 
0 473 .393 .449 
0 .629 HOC 
ee © Foo. 
0 


4Same as observed value in 1975. 


bChurchill, E., McConville, J.T., Laubach, L., and 


White, R.M., Anthropometry of U.S. Army Aviators-1970, U.S. 
Army Natick Laboratories Technical Report 72-52-CE, December 


Mogi, p. 278. 
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Means, 


Variable 


Sitting Height 
mean 
std.dev. 
range 

Functional Reach 
mean 
std.dev. 
range 

Buttock-Knee Length 
mean 
std.dev. 
range 

Leg Length 
mean 
std.dev. 


range 


Table V 


DATA SET COMPARISON 


1964 


20. 


oe 


SSI 


Zee 


24. 


ZOE 


43. 


Se 


28 


.29 


19-41 .62 


ol 


42 


26-36.31 


O09 


. 00 


73-27.81 


82 


moO 


97-49 .387 


33 


Standard Deviations, 


1969 


36. 


32 


Sal 


26. 


24. 


20. 


43. 


So), 


76 


Eo 


.43-40.61 


24 


.62 


02-35.73 


o4 


. 26 


67-28.31 


86 


nos 


15-50.02 


and Ranges 


USS 


36. 


S2e 


30. 


ou 


24. 


ZOr 


42. 


il 


61 


miles 


90-40 .00 


89 


SG 


10-34.60 


QO 


ae) ii 


70-26. 80 


og 


oo 


10-48.20 





Table VI 


DATA SET COMPARISON 








Correlations 
1964 
Sitting Functional Buttock- Leg 
Height Reach Knee Length Length 
Sitting Height 120 oO Oe ~45 
Functional Reach 1.0 -DOG 639 
Buttock-Knee Length 0 86 
Leg Length © 
1969 
Sitting Height 1.0 44 .40 245 
Functional Reach Lao G'S .639 
Buttock-Knee Length 0 86 
Leg Length .O 
1975 
Sitting Height 1720 44 .40 245 
Functional Reach © Oe: Onl 
Buttock-Knee Length 1.0 ao 
Leg Length .O 





When the assumed values described above were included 
in the data analysis, the requirements for CAPE input were 
complete. The inputs to CAPE are listed in Tables I 
through IV under Observed Data. One thousand subjects were 
generated for the 1964 and for the 1969 data sets. Together 
with the 1975 data in its original, as observed form, the 
result was three similar data sets with sample sizes of 
1000, 1000, and 936. The validity of this model's perfor- 
mance has been demonstrated with a correlation between 
empirical and generated results of 0.997 to 0.999 CBittner, 
1974]. For the 1964 and 1969 data sets, the CAPE generated 
output was analyzed to insure that the means, standard 
deviations, and correlations of the generated data were 
actually those sought. The values used as inputs to CAPE 
and those resulting from analysis of the generated data are 
Shown in Tables I through IV. Tables V and VI compare the 


summary statistics for all three data sets. 
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IV. ANALYSIS OF DATA 


As a result of the data generating procedure, all three 
data sets included a value for each of the four relevant 
measurements for each subject. The samples were then 
examined more closely to determine which subjects in each 
sample would be excluded, either by design constraints (Sth 
to 95th percentiles or 3rd to 98th percentiles) or by air- 


craft constraints specified by NAVAIRINST 3710.9. 


A. PERCENTAGES EXCLUDED 

Fortran program SCREEN was used to count the subjects 
who were too large or too small to satisfy anthropometric 
requirements. For each of the three samples, identical 
screening operations were conducted. Each subject was 
examined and excluded if any one or more of the four 
measurements under consideration was outside prescribed 
limits. Each subject's anthropometric features were tested 
against each cutoff limit. Thus an individual might be 
excluded on both functional height and sitting height but 
not on buttock-knee length and leg length. For each elimin- 
ation encountered, the subject was counted as being 
eliminated and each variable causing an elimination was 
tallied. At the end of each sample, the resultant totals 


were converted to percentages of the total sample, yielding 
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a percentage eliminated by each of the four variables and 
meeoetal exclusion resulting from the net effect of the four 
variables' eliminations. 

The prescribed limits mentioned above were set equal to 
the percentile values in the design specification (Sth and 
95th or 3rd and 98th percentiles). Then the limits were 
set, successively, to the anthropometric limits of each 
aircraft for which anthropometric coding had been estab- 
lished CNAVAIRINST 3710.9]. 

1. Design Constraints 

To screen the samples by deSign specification, the 
subjects were screened to eliminate those with any of the 
four measurements less than the Sth percentile or above the 
95th percentile of the 1964 data, the design standard [MIL- 
omp-1472BJ. Then the screening cutoff points were changed 
to correspond to the 3rd and 98th percentiles [CNAVAIR SD-24K, 
1973), and the screening process was repeated. 

It 1S interesting to note at this point that percen- 
tile values for leg length of the 1964 data do not exist 
Since that measurement was not taken in the 1964 survey. 
Still, the Navy considers leg length to be a critical dimen- 
Sion in that it is one of the measurements included in the 
Anthropometric Compatibility Assignment Program. Critical 
dimensions are also the basis of the design specification 
Since '...design limits shall be based upon a range from the 


Oth percentile to the 95th percentile values for critical 
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body dimensions' E~MIL-STD-1472B]. Since the percentile 
values were not available, values used for leg length 
percentiles were determined from the CAPE generated sample 
of 1000 subjects, described in Chapter III, DATA PREPARATION. 
The design specification percentile screens resulted 
in the maximum exclusion percentage to be expected if air- 
craft were designed to accommodate only the specified 
percentiles, and no more, for each dimension. This result 
amounts to a worst case exclusion and could occur only if 
a design exactly met the bare minimum accommodation in each 
variable, a very unlikely event, and one which would still 
not meet the requirement to accommodate 90 percent of the 
user population as specified by MIL-STD-1472B. The percen- 
tile screen results, shown in Table VII, show that cumulative 
exclusion can greatly exceed the maximum single variable 
elimination. For example, the 1969 data screened on the 
5th to 95th percentiles shows a total exclusion of 40.4 
percent but the highest single variable elimination is the 


21 percent excluded by buttock-knee length. 
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Table VII 


Design Specification Exclusions for 
the 5-95 Percentiles and the 3-98 Percentiles 


fariable 1964 Too JLo 


5th to 95th Percentile 


Percentage Excluded 


Sitting Height ae Ll ey. 7.7 
Functional Reach oes 14.6 8.4 
Buttock-Knee Length 9.1 Zao L3n2 
Leg Length onS 1a Pies 
Cumulative 26.4 40.4 26.7 


3rd to 98th Percentile 


Percentage Excluded 


Sitting Height 2.0 4.7 Zed 
Functional Reach 4.8 9.4 ono 
Buttock-Knee Length o.4 LAS 6.4 
Leg Length 4.9 4.7 6.9 
Cumulative 14.3 cee 14.3 
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The difference between cumulative exclusion and the 
largest univariate exclusion figure, which could be taken 
falsely as an upper bound, is due to the non-unity corre- 
lation between variables. Although some anthropometric 
models have described all anthropometric dimensions as a 
specific multiple of height ERoozbazar, 1979] with the 
implicit assumption of intervariable correlation equal to 
1.00, those models fail to recognize the variable propor- 
tions of people and would lead to the acceptance of the 
largest univariate exclusion as the net cumulative exclu- 
Sion. The assumption of a correlation of 1.00 amounts to 
treating the population as if the man with 95th percentile 
height must also have 95th percentile leg length, 95th 
percentile shoulder width, and 95th percentile in all 
femesrements, i.e., ‘the 95th percentile man.‘ If that were 
the case, the percentile screens would have eliminated the 
same subjects on all variables, as well as on the cumulative 
exclusion. If, at the other extreme, intervariable corre- 
lation were zero, the single variable eliminations would be 
independent and their arithmetic sum would appear as the 
cumulative exclusion. 

In reality, some subjects with leg lengths shorter 
than the Sth percentile also have sitting heights less than 
the 5th percentile. Others have short legs but sitting 
heights nearer the mean. The evidence of variable propor- 


tions along with variable sizes is found in the cumulative 
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exclusion values. Generally, the cumulative exclusion is 
less than the sum of the single variable eliminations but is 
greater than any of the single value eliminations. Only 
rounding error could cause the cumulative exclusion to 
exceed the sum of the single variable eliminations. The 
cumulative exclusion is generally as small as the largest 
single variable when only one variable resulted in any 


elimination. 


eee Aircraft Constraints 

The specific aircraft screens were accomplished with 
the same technique as used on the design specification 
screens. The cut points for each variable were the values 
corresponding to each aircraft's anthropometric restrictions 
specified by NAVAIRINST 3710.9. Then, for each aircraft, 
the three samples were examined to eliminate subjects whose 
anthropometric measurements were large enough, or small 
enough, to be excluded by the aircraft's restrictions. The 
number of subjects eliminated by each measurement and the 
cumulative exclusions were totaled and converted to a percen- 
tage of the sample. Different aircraft with the same 
anthropometric restrictions were consolidated into eighteen 
groups and the screens were executed, resulting in the 1964, 
1969, and 1975 samples being matched to all aircraft types. 
The percentage exclusions are shown in Table VIII, which 


lists aircraft types in alphabetical order. 
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Table VIII 


pxKclusloOmeoy Alrecrart Lype 


Exclusion Percentage 
Aircraft Code Variable 1964 1969 oT 
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iaote VEL (Continued ) 


Exclusion by Aircraft Type 


Exclusion Percentage 


Aircraft Code Variable 1964 1969 1975 
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Table VIII (Continued) 


Exclusion by Aircraft Type 





Exclusion Percentage 
Aircraft Code Variable 1964 1969 LOGS 
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Table VIII (Continued) 


Exclusion by Aircraft Type 


Exclusion Percentage 
Aircraft Code Variable 1964 1969 1975 


RF-4B 89 Save. ht. ©.7 1.6 0.4 
(R) - NC Gane he. ~ - - 
9 b-knee IlInth. 0.0 0.4 0.0 
- leg length ~ - = 
cumulative 0.7 2.0 0 
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B. RESULTS 

The design criterion requiring accommodation of the 5th 
and 95th percentiles results in exclusion of over 26 percent 
of the population. That result fails to meet the goal of 
90 percent accommodation. Designing to the 3rd and 98th 
percentiles, instead of the 5th and 95th, reduced exclusion 
to 14 percent. That result still does not meet the 90 per- 
cent accommodation goal. 

No aircraft type's coding eliminated as much as ten per- 
cent of the 1964 sample. The A-7 exclusion was 8.2 percent 
and was the worst case. The accommodation by all aircraft 
of more than 90 percent would seem to indicate that the 
cockpit designs included a margin of accommodation allowance 
in addition to the minimum required. Such a safety factor 
in design is not unexpected but is still technically an 
unrequired bonus. 

The exclusion percentage was not constant, but changed 
from one sample to the next. The greatest exclusions were 
found in the 1969 data while the 1964 and 1975 cumulative 
exclusions were about the same. Although cumulative exclu- 
Sions were nearly equal, the contributing factors had 
varying proportions. Sitting height eliminations decreased 
and eliminations due to buttock-knee length and leg length 
increased. 

The 1969 data exhibited other features of incongruity 


which detracted from the credibility of the three data sets 
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Comparison of Selected Variables 
and Standard Deviations 


Variable 


Height 


Leg 
Length 


Sucre 
Height 


Buttock- 
Knee Length 


NOTE: o — mean 


Figure 3 
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taken as a whole. These effects reflected the presence of 
unidentified factors, not entirely explained by the facts 
available. While height increased steadily from 1964 to 
1975, its components did not. As seen in Figure 3, which 
shows mean values and one standard deviation on either side 
of the mean, mean height increased from 69.94 to 70.15 
inches, a rate of 0.427 inches per decade. That rate con- 
firms the USAF prediction of 0.418 inch per decade [ Roebuck, 
et. al.J. However, the increased height does not appear as 
an increase in the other long bone measurements which are 
components of height. From 1969 to 1975, both mean sitting 
height and mean leg length appear to have decreased. The 
increase in height should be traceable to an increase in 
sitting height, an increase in leg length, or both, but 
neither was noted. Explanations other than erroneous data 
are possible, but not supported by knowledge of the popula- 
tion. It could be hypothesized that an increasing 
concentration of fat in the buttocks could increase both 
sitting height and leg length measurements while not 
affecting height. This effect has not been observed. 
Another possible explanation, operator induced loss of 
precision, would be accompanied by increases in the variation. 
An examination of Table V and Figure 3 indicate that standard 
deviations remained relatively constant, indicating the 


absence of excessive random error by equipment operators. 
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The inclusion of a constant error in any of the measure- 
ments in one or more data sets could account for the 
apparent inconsistency. This kind of error should be 
precluded by the use of standard procedures and calibrated 
equipment. The 1969 data were collected from measurements 
taken with a device for which no calibration was provided. 
There was no calibration procedure or equipment associated 
with the integrated anthropometric measuring device until 
after the introduction of the improved device in 1971. 
Although a means of calibrating the improved anthropometric 
measuring device was provided after 1971, the calibration 
equipment is no longer available and the date of the last 
calibration is unknown. 

The largest contributor to the apparent error is the 
diminished leg length in 1975. The reduction of nearly an 
inch is difficult to believe. Leg length is also a diffi- 
cult measurement to take accurately EMoroney, nOSOL. “ihre 
position of the subject being measured (sitting erect with 
his right leg maximally extended horizontally in front of 
him) is difficult for some individuals to attain. The 
inflexibility of some muscular legs prevents the leg from 
being straightened completely. Methods aimed at straighten- 
ing the leg can result in the hips rolling forward, causing 
the buttocks to move outward from the seat back. The 


measured leg length is then either too long (with the leg 
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straight and buttocks forward) or too short (with the leg 
slightly bent). 

The 1975 data also appeared to have been influenced by 
rounding. The frequency of whole inch and half inch measure- 
ments is greater than expected. The bias incurred in the 
mental rounding process is indeterminate and no effort was 
made to correct for it. The existence of the effect indi- 
cates a need for a procedural modification. 

The exclusion determined here were based on the consider- 
ation of only four variables. Although, of the available 
measurements, they were considered to be the most important 
with respect to accidents, injuries, and aircraft controlla- 
bility, EMoroney, 1980] they are not the only measurements 
relevant to accommodation. Administrative limitations set 
by NMPC established that no more than four variables would 
be used. The practice of collecting anthropometric data at 
the flight training entry point and the specific measurements 
taken there determined which measurements could be used. 

The inclusion of additional measurements, such as shoulder 
width, could only increase exclusions. The substitutions of 
other measurements could affect the exclusion results 


samilarly. 
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V. SUMMARY AND CONCLUSIONS 


The Navy requires airplane cockpit designs to accom- 
modate ninety percent of the Naval aviation population. In 
Support of that goal, designs are required to accommodate 
the Sth to 95th percentiles (or the 3rd to 98th percentiles) 
of critical body dimensions. The success or failure of 
aircraft designers in accommodating the population had not 
been quantified, heretofore. By examining anthropometric 
faeaetrom three populations (1964, 1969, and 1975), two of 
which were generated by the Computerized Accommodated Percen- 
tage Evaluation Model (CAPE), exclusions allowed by the 
design specification and exclusions associated with each 
aircraft type were determined. The numerical values corres- 
ponding to the percentiles named in the specification (5 to 
95 or 3 to 98) were compared to the measurements of each 
subject in each population to determine the proportion of 
subjects of each population having at least one measurement 
outside the limits of accommodation and therefor excluded. 
The results were the exclusions expected from designing to 
the specified percentiles. The aircraft exclusions were 
found by comparing aircraft anthropometric exclusion codes, 
assigned by NAVAIR, and the measurements corresponding to 
those codes, to the three populations. That process 
resulted in an exclusion percentage for each aircraft and 


each population. 
Sil 





It was found that: 
1. The design specification requiring accommodation of the 
oth to 95th percentiles excluded 26 percent of the specified 
population. 
2. Designing to the 3rd to 98th percentiles excluded 14 
percent of the population. 
3. Exclusions, both those allowed by the specification and 
those peculiar to each aircraft, varied over time. 
4. All aircraft accommodate more than the required ninety 
percent of the most recent population sample. 
5S. The proportion of exclusion caused by a given variable 
changed over time. 


6. There are inconsistencies in the data. 


A. DESIGN SPECIFICATION 

The design specification should result in at least ninety 
percent of the user population being accommodated by each 
aircraft. Reference to the 5th and 95th percentiles and to 
the 3rd and 98th percentiles of critical body dimensions 
serves only to fog the issue. The imaginary ‘95th percentile 
man,' consisting of a 95th percentile measurement for all 
variables, is assumed from the requirement to design to his 
dimensions. The tacit assumption that anthropometric 
features can be assembled and generalized to form a 95th 
percentile man is false and misleading since it assumes a 
correlation of 1.00 between anthropometric features. The 


specification would be more effective if it were to require 
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accommodation of ninety percent of the population and remain 
Silent on the accommodation of singular dimensions. Deletion 
of all reference to percentiles in the specification would be 
a productive change. However, anthropometric data should 
still be reported in percentiles for use in univariate 
applications. 

The Navy is currently enjoying the benefit of an unre- 
quired bonus in the area of user accommodation. Ninety 
percent accommodation was required, more than that is 
delivered, at least for the four variables examined, based 
on the 1975 data. Although there appears to be nothing to 
indicate that this bonus will not continue, there is like- 
wise nothing to guaranty that it shall. If a future aircraft 
design were successful in accommodating exactly the required 
ninety percent, the consequences would amount to exclusions 
greater than those being experienced on any current Navy 
aircraft and the design specification would still have been 
Satisfied. Since there is feeling in the Navy that current 
user accommodation is insufficient and since the ninety 
percent accommodation requirement is being met, either the 
Navy must change the ninety percent requirement or adjust its 
attitude toward current accommodations. The answer may be 
found in partitioning the user population into groups, each 
Seeevnieh currently flies in a particular type aircraft. 

Then specify the 'user population' to be one of those groups, 


rather than all flying personnel in the Navy, and require 
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accommodation of the entire group. This solution, while 


easy to conceive, would be difficult to administer. 


Bec OCKPIT RESTRICTIONS 

The exclusions enumerated in Table VIII indicate that 
the Navy has a minimal problem, on the global level, with 
incompatible planes and personnel. The attention that the 
problem has received would hardly seem justified by the 
miniscule exclusions indicated by the analysis. The 
expense in terms of money, man-hours, and effort reflected 
in the production of instructions by OPNAV, BUMED, CNATRA, 
and NAVAIR, together with the resultant effort necessary to 
comply with their requirements must have been justified by 
a judgement that a problem existed but the numbers in Table 
VIII do not reflect the magnitude of the problem. A reason 
that more impressive numbers don't appear in the table 
might be in the aircraft codes' inaccuracy, or in their 
failure to account for influential variables. 

The accuracy of the aircraft codes is dependent ona 
judgement factor on the part of the people who measure air- 
craft cockpits. The anthropometric codes assigned to 
aircraft are based on cockpit measurements adjusted for 
clothing and equipment worn by the user. Those adjustments 
cannot allow for all variations. The codes assigned to 
aircraft can also be affected by forecasts of the consequences 
of the codes assigned. The prospect of excluding an unsuit- 


ably large portion of the user population could have 
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influenced aircraft exclusion code assignments in the direc- 
tion of minimal effect. 

The high cost of error does not appear in the table but 
undoubtedly accounts for the attention given to such unim- 
pressive exclusion figures. The cost of losing an airplane 
Or an aircrew is considerable. That cost makes any erroneous 
assignment of personnel potentially disastrous and may 
account for the emphasis in the area of aircrew accommoda- 
non. 

Additionally, anthropometric restrictions can affect an 
officer's career rather markedly by denying him assignment 
beomune type aircraft he is motivated to fly. That fact 
makes the globally small exclusion problem a crisis on the 
individual level and explains why the attention level 


exceeds that expected from examining the exclusion figures. 


C. DATA INCONSISTENCIES 

The data indicate that the values associated with anthro- 
pometric features shifted over time. Although the exclusions 
by design specification in 1964 and those in 1975 were about 
the same (see Table VIII), the individual variables respon- 
Sible for the excluSions were of varying influence. [In 1964, 
sitting height caused more exclusions than it did in 1975, 
while buttock-kKnee length and leg length caused more in 1975 
than in 1964. The 1969 data show higher exclusion than 


either of the other two data sets. Sitting height, buttock- 
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knee length, and leg length all peaked in 1969, accounting 
for the corresponding peak in exclusions. 

The exclusion results varied erratically over time. 
There are inconsistent effects in the data and the lack of 
continuity cannot be explained with the information avail- 
able. The leg length measurement is the most likely source 
of the inconsistent effects and the discontinuance of leg 
length as a critical measurement might be beneficial if 
another measurement could be substituted. Buttock-knee 
length might be used as a substitute for leg length since 
the two variables are well correlated (r=.8). Alternatively, 
a measurement taken ina sitting position approximating a 
cockpit seating position could be substituted. Measuring 
from the heel to the seat back via the front edge of the 
seat would accomplish, more directly, the purpose presently 
being attempted by the leg length measurement. The problems 
encountered measuring leg length would be eliminated and the 
measurements of personnel and of aircraft cockpits would be 
more comparable and, consequently, more effective in matching 
personnel to airplanes. Such a measurement of ‘functional 
leg throw' would be a Special purpose dimension not included 
in previous surveys. Previous surveys, however, were not 
faced with the special problem now at hand and the cost of 
designing and implementing this new measurement could be 


recovered in the prevention of a single aircraft accident. 
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D. DATA LIMITATIONS 

Of all the anthropometric dimensions which might be con- 
sidered, the Navy collects data on eight anthropometric 
features (including weight) considered relevant to flight 
safety. Of those eight variables, only four are considered 
in the process of assigning personnel to aircraft. The 
limitation to four variables is not necessarily optimum with 
respect to flight safety since the administratively driven 
maximum of four variables was the effective constraint. The 
inclusion of additional variables, if they were found to be 
relevant, could only decrease accommodation further. The 
decrease could be dramatic if the additional variables were 
poorly correlated to those already considered. 

The problems arising from variables in combination (but 
not visible in single variable evaluation) have been recog- 
nized [CNATRAINST 13520.1, 1980] but no description of the 
interaction between variables is included in the Anthropo- 
metric Compatibility Assignment Program. For example, the 
minimum functional reach and leg length necessary to reach 
controls varies as the position of the seat changes. The 
position of the seat is determined largely by the sitting 
height of its occupant. An individual with a short sitting 
height must raise the seat high enough to allow vision over 
the nose of the aircraft. As the seat rises, it also moves 
aft slightly, causing the rudder pedals to be adjusted aft 


and also requiring a longer functional reach to enable the 
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man to reach the controls. Thus, sitting height determines 
in part the effective minima for leg length and functional 
reach values, making required leg length and functional 

reach distances functions of sitting height. The values used 
for anthropometric coding, however, are determined irres- 
pective of seat position changes. An additional 
anthropometric code should be established to identify those 
personnel who meet all the limits but whose combination of 
anthropometric features could create problems. 

The procedures used to take measurements in the 1969 and 
1975 data sets were deficient in some respects. Most notably, 
the failure to calibrate the integrated anthropometric 
measuring device was unfortunate. The apparent rounding of 
measurements also induced error. The calibration deficiency 
can be corrected by instituting a requirement to document 
periodic calibration and adjustments, as is required by the 
Preventive Maintenance System. The rounding error could be 
eliminated by converting the measuring device to provide 
digital displays of measurements. Better, but also more 
expenSive, an automated measuring device could record 
measurements without requiring the operator to read and 
record numbers. Such a device would also prevent the oper- 
ator from rounding measurements or inserting any bias into 
the data. 

The necessity to draw on outside sources for data con- 


cerning leg length, a critical body dimension about which 
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there is no data in the design specification (1964 data), 
points out a basic deficiency in the specification. 
Designers do not consider leg length in the design of 
cockpits. Instead, they use buttock-kKnee length and pop- 
Picateheight sitting, which are components of leg length. 
The assignment process would be more consistent if similar 


variables were used for exclusion codes. 


E. CONCLUSIONS 

The data collection system which supports the Anthropo- 
metric Compatibility Assignment Program consists of 
measurements taken on aircraft cockpits as well as the 
equipment and procedures used to measure aviation personnel. 
Many inconsistencies in personnel measurements have been 
noted even though those measurements are an essential ingre- 
dient in the successful operation of the program. The 
program will operate no more effectively than its weakest, 
or limiting, factor. That limiting factor appears to be the 
measurement of personnel. Inconsistencies in the data appear 
to be results of measurement bias. The appearance of larger 
personnel in 1969 and corresponding higher exclusions in 
that population are more likely the changes of measurement 
bias than any real effect in the population. 

Measurement error could be reduced by: 
1. Periodically calibrating the measurement equipment. 


2. Randomly remeasuring personnel and comparing the original 
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and repeated measurement. If significant discrepancies are 
noted, corrective action should be taken. 

3. Replacing the integrated anthropometric measuring device 
by one using digital measurement displays, or automated 
measurement recording. 

The consequences of failing to improve the Navy's 
personnel anthropometric measuring procedure will amount to 
negating the Anthropometric Compatibility Assignment Program 
and accepting the aircraft accident rate, and ejection injury 
rate which were sufficiently undesirable to lead to the 
program's establishment. Additionally, erroneous measurement 
can adversely impact on officers' career development patterns 


and consequently on officer personnel retention. 
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APPENDIX A 


Definitions of Measurements 


Buttock-Heel Length -see Leg Length 


Buttock-Knee Length -the horizontal distance from 
the rearmost surface of the buttocks to the front of 
the kneecaps, with the subject sitting erect. 


Buttock-Heel Length -see Leg Length 


Functional Reach -the distance from the wall to 
the tip of the thumb measured with the subject's 
shoulders against the wall, his arm extended forward, 
and his index finger touching the tip of his thumb. 


Height -vertical distance from the 
floor to the top of the head with the subject standing 
erect in bare feet. 


Knee Height Sitting -the height from the footrest 
surface to the musculature just above the knee. 


Leg Length -the distance from the base 
of the heel to a wall against which the subject sits 
erect with his leg maximally extended forward along 
the sitting surface. 


Sitting Height -the height, from the sitting 
surface, to the top of the head, with the subject 
Sitting erect. 


Stature -see height 


Popliteal Height Sitting -the height of the underside of 
the upper leg above the footrest surface. 
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APPENDIX B 


Anthropometric Codes as Percentiles 
of the 1964 Population 


Variable Code Measurement Equivalent 
Interval Percentile 
SITTING HEIGHT 9 40.0-41.0 99+ 
8 39,5-39.9 99+ 
V 39.0-39.4 98-99 
6 38 .5-38.9 96-97 
5 38.0-38.4 91-95 
4 35.0-37.9 15-90 
3 34.0-34.9 3214 
2 30 0=30 09 £3 
1 32.5-32.9 <1 
0 92.0232 .4 <1 
FUNCTIONAL REACH 9 $27.9 <1 
8 28.0-28.4 <1 
fk 28.5-28.9 La 
6 29.0-29.4 3-6 
5 205 = 304 fais 
4 30.5-30.9 26-39 
3 Se Oa ees 40-50 
p 31.5-32.4 Sao 
il 32.5-33.9 76-94 
0 SO >95 
BUTTOCK-KNEE LENGTH 9 > 28.0 >99 
8 2770223 ..0 >99 
ii 26.5-26.9 >98 
6 26.0-26.4 97-98 
5 25.5-25.9 91-96 
4 25.0-25.4 81-90 
3 24.0-24.9 50-80 
2 23 0=25°.9 15-49 
i 22 .0-22.9 Day 
0 S$ 21.9 <2 
LEG LENGTH 9 49.0-50.0 >99 
8 48.0-48.9 98-99 
i AT MO=47 49 95-97 
6 46.0-46.9 86-94 
5 45.0-45.9 71-85 
4 ASeO=44..9 34-70 
3 40.0-42.9 Bead 
2 39,0-39.9 ieee? 
1 38.0-38.9 <1 
0 36.0-37.9 <1 








Appendix C 


Program Listing for CAPE Model 
Used to Generate the 1964 and 1969 Data 








PROGRAM GENSAMP (FMEAS eUME AS «SAMPLE oe INFUT SOUTRUT.« 


] TAPF 1=FMEAS« TAPE 2=UMEAS« TAPE J=SAMPLE 0 
2 TAPE S=INPUT + TAPFA=QUTPUT) 
Cc 
Cc eee MODULE A = MONTE CARLO SUMULATTON 90° 
GC #27 MCSM READS IN FACTOR MATRIX UF REQUIRED ANTHROPIMETRIC 
c eee MEASUREMENTS FROM FITHER DISK FILES OR INTERACTIVE 
Cc #ee@ TFOMINAL. IT ALSO ALLOWS USFRS TO PRINT» EDIT OND SAVE 
Cc #@0? FILES. THE PROGRAM THEN COMPUTE RANDOM MULTIVARIATE NORMAL 
Cc eee FFATURE VECTORS, AT THE ENJ OF THF PROGPAM IT CALLS 
c @@@ SURPOUTINE TRANS TQ PERFORM NECESSARY DATA TRANSFORMATION 
Cc @@2 OF ANTHROPOMETRIC SAMPLES AND STORFS THE RESULTANT 
Cc @@@ MEASUREMENTS ON DISK FILES FOR CFEWSTATION EVALUATION 
Cc eee MONULE, 
Cc ‘ 
COMMON/SMCSY/CORR (78) eANTM( 2012) 0 ¥(25) oNMEAS oKNe 
1 SQROOT (78) eX (25) eIUeNTMEAS es ICOUNT® 
2 INAME (4) eI TMEAS -NIMEAS eNTEMP 
COMMON/MISC/IOINs LOOUTeICARD 
Cc 
Cc NMEAS = NUMBER OF ANTHROPOMETRIC MEASUREMENTSe CURRENTLY 12 
c NSAMP - NUMBER OF OPERATOR SAMPLES 
C DATA FILE DESCRIPTION? 
e NAME PROGRAM REF, TAPE# FO/UN CONTENTS 
Cc FMFAS -oc-- l FOX ANTM (2 eNMEAS) oe CORR (KN) 
Cc UME AS NTME AS 2 UN NMEAS eKNe ANT“ 0 CORR 
Cc SAMPLE NTE MP 3 UN INAME (4) oVoXk 
Cc INPUT crore S -- were or ees eccereece 
C OUTPUT eoece 6 -- eee oo cess cee eeres 
Cc 
DATA ININeOOUT sNUCREWeNTEMP eNTMEASONIMEAS/S5 p62 Ge3e4e 7/ 
DATA MAXFUNC/4/ 
DATA MMFAS eNTMEAS eoNIMEAS/S120204/ 
lLeTCOUNT/O/Z 
DNATA ANTM/242®0/4CORR/7382®0/ +I NAME /S4@)] OH / 
DATA MAXOPT eMAXANTHeMAXACT/3 0693/7 
Cc 
C *#®@COMPUTFE (KN) THE LENGTH OF THE LINEAR ARRAY CORR(AKN) 
Cc 
1OOUT=0 
ISAMP=0 


C *@@MCSM MAIN MENU 
100 CONTINUE 
IOIN=5 
WRITE (60290) 
CALL INPT (lel eNVALSZeNMEAS eo XVAL eC IEOF) 
GO TO (13041000110 +12062200120)s EOF) 
110 CONTINUE 
CALL HFLP (1) 
GO To 100 
120 CONTINUE 
WRITE (&+340) 
60 TO 100 
130 CONTINUE 
IF (NMFASSLT.1 -OReNMEAS.GT.MMFAS) GO TO 120 
KN=NMF AS® (YMEAS01)/2 
TINPT=0 
Cc 
C ##0@ GENERATE OPERATOR SAMPLE FROM MFANSeSTD VEV & CORK MATRIA 
Cc 
140 CONTINUF 
IOIN=S 
C *®P®OANTHROPOMFTRIC FUNCTIONS 
150 CONTINUE 
WRITE (#2300) 
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CALL INPT (leleNVALSZeIOPToeXVAL eo IEOF) 
GO TO (180+140+1604170+280+0170)+ IEOF +1 
160 CONTINUE 
CALL HELP (2) 
GO TO 150 
170 CONTINUE 
WRPITF (62340) 
GO TO 150 
140 CONTINUE 
IF (TOPT.LT.~1-ORe-IOPT.GTeMAXANTH) GO TO 170 
C INPUTSEDITsSAVE ePRINTeGENERATE OPER SAME+END 
G60 TO (1904200+220924042605280) > IOPT 
190 CONTINUE 
C ##@8 INPUT ANTHROPOMETRIC MEASUREMENTS 
CALL ACCPT (ISAMP) 
LINPT=] 
G0 TO 140 ‘ 
200 CONTINUE 
C 8°ENIT ANTHROPOMETRIC DATA FILE 
IF (TINPTLEQ.1) GO TO 210 
WRITE (69350) 
GO TO 140 
210 CONTINUE 
Call @&nrt 
GO TO 140 
220 CONT INUF 
C 8#82S5AVE ANTHROPOMETRIC DATA FILE ON UMEAS 
IF (ITSasP.EQ.21) GO TO 230 
WEITE (&0320) 
GO TO 140 
230 CONT INUE 
CALL SAVE 
G60 TO }40 
240 CONT INUE 
C 888POQOINT ANTHROPOMETRIC DATA 
IF (TINPT.EQ.1) GO TO 250 
WRITE (42350) 


GO TN 140 
250 CONTINUE 
CALL PRINT 
GO TO 140 
c 
C ##2 GENERATE OPERATOR SAMPLE FROM ANTH MEAS DATA 
c 


260 CONTINUE 
IF (TS4MP.NE.1) GO TO 270 
CALL OPSAMP (ITSAMP) 
60 TO 140 
270 CONTINUE. 
W#PITE (459330) 
GO TO 140 
C ##e2TFRMINATION OF MCSM 
280 CONTINUE 
WOTTE (60310) 
STOP 


290 FORMAT ("" ENTER NUMBER OF MEASUREMENTS (1-12) ?-<") 

300 FORMAT (*% GENERATE FUNCTIONS(1<INPUT §2*EDIT33-SAVES4=PRINT Sy 
+ "“S-GENERATE $6-END) o-') 

310 FORMAT (* END OF OPERATOR SAMPLE") 

320 FORMAT (** UNABLE. TO SAVE MEASUREMENT DATA = FATAL ER¥ORS") 

3390 FORMAT ("* UNARLE TO GENERATE SAMPLE = FATAL ERRORS") 

340 FORMAT (" INVALID RESPONSE") 

350 FORMAT (" ANTHROPOMETRIC DATA MUST BE: INPUT BEFORE TmE EDIT '"e 
+ “FUNCTION CAN SE USEN") 

WAN FORMAT. ¢*§ ANTHROPOMETRIC NATA MIST RF TNPUT HFFORE THRE PRINT*. 
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PAAANAAANANAAARAAANAANAARAARDA 


100 


120 


130 


140 


150 


ano - 


160 


170 


+ “FUNCTION CAN 3E USED!) 
END 
SUBROUTINE INPT (ITYPE «I VALSZeNVALSZsIVAL sXVAL e+ IEOF ) 
COMMON/MISC/IOINs IOOUT e ICARD 
FREE FIELD INPUT ROUTINE -= READS IN INPUT FROM THE KEYROARD IN 
FREE FIELO INPUT 
INPUT wom 
ITYPE = SPECIFIES Tre TYPE OF DATA TU BE ENTERED 
1ZINTEGER 2=REAL 3=ALPHANUMERIC 
. IVALSZ = THE NUMBER OF ITEMS EXPECTED TO BE ENTERED 
OUTPUT om= 
NVALSZ = THE NUMBER OF ITEMS ENTERED BY THE USER 
IVAL = ARRAY CONATINING INTEGER OR 
ALPHANUMFRIC DATA SUPPLIED BY THE USER 
XVAL ~ ARRAY CONTAINING REAL DATA SUPPLIED Br THE 
USER 
IEOF = END OF FILE FLAG 
Y=IF BLANKS Or A CARRIAGE CONTROL IS THE 
ONLY INPUT 
2= IF A QUESTION MaRK IS ENTERED 
3=IF INPUT IS INVALID 
GS=IF A $ HAS BEEN ENTERED 
DIMENSION ICHAR(80) eI VAL (IVALSZ) «XVAL(IVALSZ) 
ISTOP=9 
IF (IV4tSZ.6E.0) GO TO 100 
ISTOP=] 
IVALSZ=IAB8S(IVALSZ) 
CONTINUE 
IQUOTE=0 
NVALSZ=0 
IF OF=0 
IF (IOOUT-EQ.1) GO TO 120 
CO 110 J=1,80 
ICHAR (I) 214 
CONTINUE 
READ IN TRE INPUT 
READ (IOIN+330) ICHAR 
IF (IOIN.NE.S) WRITE (64330) ICHAR 
IF (FOF C(IOIN).NE.1) GO TO 130 
TE OF =] 
RETURN 
CONTINUE 
I1=ICARD 
GO TN 140 
CONTINUE 
INITIALIZE THE CARD COLUMN COUNTERS 
T1l=1 
CONT INUF 
12280: 
CONT INUF 
IF (ISTOPeEQ.1 -AND-IVALSZeEQeNVALSZ) RETURN 
SEFARCH FOR THE NEXT WORD 
IF (I11.LE-I2) GO TO 160 
SEARCH FOR THE FIRST NON-BLANK CHARACTER (START OF WORD) 
CHECK IF EXCESS HAS BEEN ENTEREO 
IF (NVALSZ.LE-IVALSZ) RETURN 
TEOF=5 
RETURN 
CONTINUE 
00 170 IT=I1l-eJ2 
IF (ICHAR(I) EQe1H ) GO TO 170 
IFIRST=] 
GO TO 180 
CONTINUE 
IF (NVALSZ.EQ.0) IFOF=] 
CHECK IF EXCESS HAS REEN ENTERED 
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180 


190 


200 


210 


220 


230 


240 


250 


260 


270 


280 


290 


IF (NVALSZ.GT.IVALSZ) JEOF=S 

RETURN 

CONTINUF 

SFARCH FOR LAST NON-8LANK CHAPACTER (END OF wORD) 
IF (CICHARC(CIFIRST) -EQ.1H?) IEOF=2 

IF (ICHAR(IFIRST) -EQ.1HS) IFEOF=4 

IF (IFEOFeNE.0) RETURN 

IF (ITYPE-EQ.4) GO TO 300 

IF (ITYPE-NE.3) GO TO 190 

IF (ICHARP(CIFIRST) NETH" .AND.ICHAR(IFIRST) -NE.1H*) GO TO 190 
IFIRST=IFIRSTo]) 

IQUOTE=}1 

CONTINUE 

NO 220 I=IFIRSTsI2 

IF (IOQUOTE.EQ.0) GO TO 200 

IF (ICHAR(T) EQ. 1H" .OR.ICHAR(I) .€Q.1H') GO TO 210 
GO TO 220 

CONTINUE 

IF (ICHAR(T) NEolH eANDLICHAR(IT) eNEsiHs) GO TO 220 
CONTINUE — 

ILAST=J-] 

60 TO 230 

CONTINUE 

ILAST=]? 

CONTINUE 

NVALSZ=NVALS7Z¢] 

RESET THE FIRST CHARACTER TO RE SEARCHED FOR THE NEXT WORD 
INCREMENT THE NUMSER OF WORDS FOUND 


“TIS=ILAST+? 


IF (NVALSZ.GT.IVALSZ) GO TO 150 

CHECK IF INPUT IS ALPHA=-NUMERIC 

IF (ITYPE-£Q.3) GO TO 310 

INPUT TS NUMERTC 

CHECK JF NUMSER IS SIGNED 

ISIGN=} 

IF CICHARC( IFIRST) NE elH=eANDeICHAR( IFIRST) eNEe1H+) GO TO 240 
IF (ICHARCIFIRST) eEQelH=) ISIGN=-1 
IFIPST=IFIRST] 

CONT INUF 

INUM=0 

INEC=0 

IDECF=n 

DO 260 J=IFIRSTsILAST 

IF (ICHAR(T) .NE.LH.) GO TO 250 
IDECFHINECF o} 

INEC=1 

G40 TO 2760 

CONT INUF 

CALL ICONV (ICHAR(I) eIADDoIEOF) 

IF (IEOF.EQ-3) RETURN 
INUM=IN'UM#104T A400 

IF (IDECeNE.0) IDEC=IDECe]1 

CONTINUE 

IF (IDOECFeLT.2) GO TO 280 

CONT INUF 

IFOF=3 

RETURN - 
CONTINUE 

INUM=INUM*®ISIGN 

XNUM=FLOAT (INUM) 

IF (CITYPE.EQ.1.- AND. IDEC.LE.1) Ore (I TYPE.EQ.2)) GU TO 290 
IE OF =3 

RETURN 

CONT INUF 

TE (CTTYPF.FQ.1) TVA! (NVALSZI STFA CXNUM) 
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IF (LEYPReEUlC) AVAL (NVALSZ) =XNUM 
ICARD=ILAST*) 
GO TO 150 

300 CONTINUE 
ILAST=IVALSZ#10 

310 CONTINUE 

c INPUT IS ALPHANUMERIC 
DO 320 J=leIVALSZ 
IVAL (J) 10H 

320 CONTINUE 
L=ILAST<IFIRST+]} 
NVALSZ=(L=-1])/101 
IF (LeGTelO*#IVALSZ) GO TO 270 
FNCODE (Le330eIVAL (1) ) CICHAR(IT) sI=IFIPSTsILAST) 
ICARD=1] 
RETURN 


iS 
330 FORMAT (80A1) 
END 
SUBROUTINE ICONV (JUCHAR IIe IEOF) 
Cc THIS FUNCTION CONVERTS ALPHA NUMBERS TO NUMBER NUMBEKS 
NDIMENSTON JNUM(10) 
DATA JNUM/1HON eH] eLH2eLH3s1LH4elHSelH6s]H/s1lHBe1H9/ 
DQ 100 Yy=lel0 
IF (JCHAR.NE.JNUM(J)) GO TO 100 
TI=J-)] 
RF TURN 
100 CONT INUE 
c INVALIN IN®SUT --= NOT A NUMBER 
IEOF =3 
RF TURN 
END 
SUBROUTINE ACCPT (ISAMP) 


C ##2 INPUT ANTHROPOMETRIC DATA 
C ### MEANS*e STD DEV>e AND CORRELATION MATRIX 


COMMON/MCSM/CORPR (78) sANTM(2¢12) oV (25) eNMEAS oKNe 
} SORNOT (78) 9X (25) sJUeNTMEAS eo ICOUNT 9 
2 INAME (4) s ITMEAS eNIMEASsNTEMP 
COMMON/MISC/IOINe IOOUT »I CARD 

DIMENSTON 9SUFF (3) « TEMP(12) 

DATA BUFF /3*#0/ 

ISAMP=0 

ICOUNT=90 

INLINE=0 

IOIN=S 


@e@READ MODE - 1L=INTERACTIVE 2=FILE 


KH Oan 


00 CONTINUE 
WRITE (6e610) 
CALL INPT (1e1eNOINPT+«IMODE oXINPT e+ IEOF ) 
GO TO (13059110¢120011094900110)+ IEOF*) 
C INVALID RESPONSE 
110 WRITE (6650) 


GO TO 100 

C HELP 

120 CALL HELP (4) 
GO TO 100 


130 IF (IMONE-EQ.1) GO TO 200 
IF (IMONE+«NE.2) GO TO 110 
c 
C ##@ FILE TYPES 1=FORMATTED 2=UNFORMATTED © 


Cc \ 
140  CONTINUF 
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WRITE (62620) : 
CALL INPT (IelsNOINPTeTFFORMeXINPT Ss IEOF) 
GO TO (17091509160415094900150) 5 IEOF*I 
150 WRITE (69650) 
60 TO 140 
160 CALL HELP (5) 
60 TO 140 
170 IF (IFFORM,EQ.1I) GO TO 190 
IF (IFFORM.NE.2) GO TO 150 


REwIND NTMEAS 
Se2RFAD THE NUMBER OF ANTHROPOMETRIC MEASUREMENTS» 
#OOCORRELATION COMPONENTS 
READ (NTMEAS) NMEAS oKNo INAME 
IF (ENF (NTMEAS).NE.I) GO TO 180 
WPITE (62670) 
60 TO 490 
C ®#8*7REFAD THF MEAN AND STANDARD DEVIATION VALUES 
180 RFAD (NTMEAS) (CCANTM(10+eJ0) es IO=1]1 ¢2) »JOSI sNMEAS) 
C ##2#READ THE CORRELATION MATRIX VALUES 
READ (NTMEAS) (COPR( IO) +102] 9KN) 
REWIND NTMEAS 
WPITE (62590) 
ISAMP2=] 
RF TURN 
Cc #8 SET INPUT LOGICAL UNIT TO I FOR FILE FMEAS 
190 IOIN=I 
200 CONTINUE 
c 
C *##0ACCEPT OPERATOR SAMPLE FILE NAME 
Ec 


Cc 
c 202 UNFORMATTED READ 
c 
c 


CALL NAMEFL (IERROR) 
ICOUMNT=7COUNT*¢IERROR 
INLINESINLINE®] 
Cc #2 INPUTS MEANS AND STANDARD DEVIATIONS INTO ARRAY ANT™ 
IF (CIMONE-EQ.1) WRITE (69500) 
IF CIMOMPE-NE.1]) WRITE (65510) INAME 
MAXINPT=2 
DO 320 TROW=I »NMEAS 
210 CONTINUE 
IF CIMODE-EQ.1) WRITE (69530) IROW 
220 NAMT=Q 
230 CONTINUE 
ITOTIN=MAXINPT=NAMT 
ICOL=NAMT +] 
CALL INPT (2eITOTINsNOINPT os IVAL sANTM(ICOL » IROW) » ICOF ) 
IF (IFOF.EQ.SeOR-TEOF.EQ.0) NAMT=NAMT*NOINPT 
IF (IMONE-EQ.1) GO TO 240 
INLINESINLINE*I 
IF (NAMT.LEQ.MAXINPT) GO TO 290 
IF (IEOF.EQ.-1) GO TO 480 
WRITE (62540) IROW 
IF (IEOF.EQ.5S) GO TO 270 
WRITE (A700) INLINE 
ICOUNT=I] COUNT] 
GO TO 320 
C#*% INTERACTIVE IEOF PROCESSING 
240 CONTINUE 
GO TO (280421022503 26054909270)s IEOF*! 
250 CONT INUF 
CALL HELP (6) 
GO TO 210 
260 CONTINUE 
WEITE (62650) 
GN_TO 2710 
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270 


280 


290 


300 


310 
320 


AAaNOY 


330 


340 


359 


C #92 


360 


370 


380 


390 


CONTINUE 
IOLIST=NOINPT-ITOTIN 

WOITE (6*600) IOLIST 

IF (IMODEeNE.1) ICOUNT=ICOUNT+] 
GO TO 290 

CONTINUE 

IF .(NAMT.GE.2) GO TO 290 

WRITE (69550) 

GO TO 230 

CONTINUE ; 

IF (NAMT.GT.2) NAMT=2 

IRERR=0 

DO 300 [COL=]+«NAMT 

CALL REASON (1eICOLsIROWsANTM(1eIROW) sANTM(29IROW) s IFATAL) 
IRERR=ITRERR*IFATAL 

CONT INUF 

IF (IRERR.EQ.0) GO TO 320 : 
IF (IMONE-NE.1) GO TO 310 

WRITE (69630) IROW 

GO TO 220 

ICOUNT=ICOUNT*IRERR 

CONT INUE 


242 ANTHROPOMETRIC CORRELATION MATRIX VALUES ARE aCCeErTED 
#¢% AND PROCESSED. 


IF (IMOMNE.EO.1) WRITE (62560) © 

90 460 [ROw=] sNMEAS 

CONTINUE 

NAMT=0 

MAXINPT=NMEAS+1] =I ROW 

CONTINUE 

IF (IMONE-EQ.1-AND.IROW.EQ.NMEAS) wRITE (69580) TRO eMAXINPT »NME4S 
IF (IMODE-EQ.1-AND.IROW.WNESNMEAS) WRITE (69570) IXOWeMAXINPT es IROWs 


INMEAS 


CONT INUF 

ITOTIN=MAXINPT=<NAMT 

ISTART=IROW+NAMT 

CALL INPT (2eITOTINeNOINPTeIVALeTEMP(ISTART) sIEOF) 
IF (IEOF EQ.-0-OR-IEOF2EQ.5) NAMT=NAMT*NOINPT 
IF (IMONE.~£Q.1) 60 TO 360 

INLINFESINLINEo] 

IF (NAMT.EQeMAXINPT) GO TO 410 

IF (IEOF.EQ.90) GO TO 350 

IF (IEOF E921) WRITE (62710) 

IF (IEOF.EQ21) GO TO 470 

IF (IENF-EQ.S) GO TO 390 

ICOUNT=I COUNT} 

WRITE (69700) INLINE 

GO TO 4A0 

INTERACTIVE INPUT ERROR PROCESSING 

CONTINUE 

GO TO (4004380+370238054902390)5 ILEOF?)} 
CONTINUE 

CALL HELP (7) e 
GO TN 340 

CONTINUE 

WRITE (62650) 

GO TO 340 

CONT INUE 

IOLIST=NOINPT-ITOTIN 

WRITE (69600) IOLIST 

IF (IMONDE e«NE.1) ICOUNT=ICOUNT?]1 

NAMT=MAXINPT 

GO TN 410 


Og 





400 


410 


420 


430 


440 


450 
440 
470 


480 


490 


500 
sic 
520 
530 
540 
550 
S60 


570 
540 
590 
600 
AlN 
§20 
630 
640 
650 
660 
670 


640 
690 
700 
710 


CONTINUE 


IF (NOTNPT.EQ.ITOTIN) 


GO TO 410 


ITOTIN=MAXINPT=NAMT 


WPITE 


(69590) 


ITOTINe I20W 


GO TN 350 
CONT INUF. 
ILAST=]ROW+NAMT=-1 


IRERR=0 


DO 420 ICOL=IROW-+ILAST 


C4LL RFASON 


(2esITCOLeIROWs TEMP (ICOL) sOeJTFATAL) 


TRERR=TRERR+IT FATAL 


CONT INUF. 

IF (IRERR.EQ.0) GO TO 440 

IF (IMONE-EQ.2) GO TO 430 

WRITE (69640) IROW 

GO TO 310 

ICOUNT=ICOUNT¢IRERR 

GO TO 440 : 
CONTINUE 


NO 450 TCOL=IROWsNMEAS 
N=NDX (TROWsITCOL) 
CORR(N) =TEMP(TCOL) 
CONTINUE 

CONTINUE 

CONTINUE 


IF (ICOUNT.EQ.0) 


ISAMP=]} 


IF (IMONE+NE.1) WRITE (69680) ICOUNT 
GO TO 490 


WRITE 
RETURN 


(69660) 


CONT INUF 


RE TURN 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 


1E UPPER "7 TRIANGULAR FORM."/" NOTES 


(" ENTER PAIRS OF MEANS AND STANDARD DEVIATIONS=<="") 
(/7** DIAGNOSTICS OF FORMATTED RKESD "7" OPER DESC: "e4Al10/) 
(77" CORRELATION/SANTHROPOMETRIC MEASUREMENTS?" ) 

( MEASUREMENT = "sI2e't Pema) 

(" MEASUREMENT -—"5]2) 

(" ENTER STANDARD DEVIATION =") 

(" ENTER ANTHROPOMETRIC CO*¥RELATION MATRIX DATA IN A ROWWIS 
ALL DIAGONAL ELEMENTS MUST" s 


" FOQOUA4L 1.000") 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 


(SH POW sI2o"("sT2e" ELEMENTSe COLUMNS "olTest="sTee%)") 
(SH POW sI2e"(",T2e" ELEMENTs COLUMN "oTes'')"*) 
(2lq ENTER THE REMAINING sI2s1#H ELEMENTS FOR ROW oT2) 


(* THE LAST "elI2s'* INPUTS HAVE REEN IGNORED") 
(" INPUT MONE (1=INTERACTIVES 2=F ILE) -—") 

( READ MODE (1=FORMATTED’ 2=UNF ORMATTED) -='") 

(" RE-ENTER DATA FOR "9 2) 

(" RE-ENTER DATA FOR ROW "sJ 3) 

(" INVALID RESPONSE") 

(" UNEXPECTED EOF ON FMEAS FOR MEAS 4] 2) 

(" FILE UMEAS IS NOT ATTACHED TO THE PROGRAM'"*/ 


" UNFORMATTED READ ON THIS FILE IS ABORTED -r="9//) 


FORMAT 
FORMAT 


FORMAT 
FORMAT 


END 


(* TOTAL ERRORS ON INPUT VATA "4]5) 

(" MEANsSTD DFV & COPRELATION DATA READ") 
(* JINVALIN DATA AT LINE ""o15) 

("* UNEXPECTED ENO OF FILE") 


SUbBPOUTINE EDIT 
COMMON/MCSM/CORR (78) sANTM (2912) 9 V(25) oaNMEASSKNo 
SQRONT (78) eX (25) eIUeNTMEASs ICOUNT » 
INSME (4) gpI TMEASsNIMEASsSNTEMP 
COMMON/MISC/IOINe IOOUTs I CARD 
DIMENSION SUFF (3) 
Asta BIFEZ ISOS 
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oe ke ke) 


100 


120 


130 


KOOaAND 


159 


160 


170 


160 


190 


200 


230 


NAD 


#2? MONIFICATIONS ON MEANSsSTANDA4RD DEVIATIONS» 
#@% AND ANTHROPOMETRIC CORRELATIONS 


IOIN=5 

CONTINUE 

WRITE (69350) 

CALL INPT (leleNVALSZsIOPTsXVALsIEOF) 
GO TO (130910091109120+340+120)s IEOF*) 
CONTINUE 

CALL HELP (8) 

GO TO 100 

CONT INUE 

WRITF (69360) 

60 TO 100 

CONT INUF 

IF (ITOPT.LT.1.0OR.~IOPT.GT.5) GO TO 120 
GO TO (1409140%240+330+340)s IOPT 


##% TCOL=1 FOR MEAN PROCESSING 
#2 TCOL=2 FC STANDARD DEVIATION PROCESSING 


CONTINUE 

ICOL=IO0PT 

CONT INUE 

MAXINPT=2 

NAMT=0 

CONT INUF 

IF (ICOIL.EQ.1) WRITE (62380) 

IF (ICOL.E9.2) WRITE (62390) 
CONTINUF 

IVALSZ=MAXINPT-NAMT 

J=NAMT +s] 

CALL INPT (2seIVALSZeNVALSZeIVAL 9 SUFF (J) s IEOF ) 
GO TO (220+180+200019091009210)>s IEOF?) 
CONTINUE 

WPITE (60370) 

GO TO 340 

CONTINUE 

WPITE (6+*360) 

GO TO 160 

CONTINUF 

IF (ICOL-eEQ.1) CALL HELP (11) 

IF (ICOLe~EQe2) CALL HELP (12) 

GO TO 140 

CONTINUE 

IOLISTSENVALSZ-IVALSZ 

WEITE (68400) IOLIST 

CONTINUF 

IROW=SIIFF (1) 

NAMT=NVALSZ*NAMT 

IF (NAMT.GEeMAXINPT) GO TO 230 

IF (ICOL.EQ.1) WRITE (6420) JROW 
IF (ICOIL-EQ.2) WRITE (69430) JROW 
GO TON 170 

CONTINUE 

CALL REASON (LeICOL+IROWs BUFF (2) 905 ITF ATAL) 
IF (IFATAL.EQ.1) GO TO 150 
ANTM(ICOL»s ITROW) =BUFF (2) 

Go TO 100 


#02 PROCESSING FOR ANTHROPOMETRIC COFRELATIONS 
CONTINUE 


MAXINPT=3 
NAMT=0 
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250 


260 


270 


220 


290 


300 


310 


320 


c 
e 


WOaAAaAA 


340 


350 
340 
370 
380 
390 
400 
410 
420 
430 
440 
450 


CONTINUE 

WPITE (452410) 

CONTINUE 

IVALSZ=MAXINPT+NAMT 

J=NAMT +] 

CALL INPT (2sIVALSZeNVALSZeIVAL e BUFF (J) - LEOF) 
GO TO (310+270+290+280+100%300)s IEOF*) 
CONT INUE 

WPITE (62370) 

GO TO 250 

CONT INUE 

WRITE (69360) 

GO TO 250 

CONTINUE 

CALL HELP (13) 

GO TO 250 

CONTINUE 

IOLIST=NVALSZ=-IVALSZ 

WPITFE (62400) JIOLIST 

GO TO 320 

CONTINUE 

ICOL=BUFF (]) 

IPOW=BUFF (2) 

NAMT=NAMT *NVALSZ 

IF (NAMT.GESMAXINPT) GO TO 320 

IF (NAMT.CEQe]1) WRITE (62440) ICOL 

IF (NAMT.EQ.2) WRITE (62450) ICOLsTIROW 
GO TO 2760 

CONTINUE 

CALL REASON (3+I1COLs IROW+BUFF (3) sOeTFATAL) 
IF (IFATAL.LEQe1) GO TO 240 

CALL RFASON (2eI1COL s+ IPOWseBUFF (3) sUsIFATAL) 
IF (CIFATAL.LEQ.1) GO TO 240 


#001 0AD THE VALUE OF BUFF (3) INTO T7E APPROPRIATE 
#@ePOSITION WITHIN THE CORRI(N) ARRAY 


N=NOX (T20wWe ICOL) 
CORR (N) =AUFF (3) 
GO TO 100 


#@% PRINT OPTION 


CONT INUE 

CALL PRINT 

GO TO 100 

#e2 END OF EDIT MODE 
CONTINUE 

RETURN 


FORMAT ("' EDIT MODE (1-MEANS32-STU3 3-ANTH.CORR$ G-PRINTS5<-END) -='") 
FORMAT (" INVALION RESPONSE") 
FORMAT (' NO INFORMATION ENCOUNTERED") 
FORMAT (* INPUT INDEX AND MEAN VALUE") 
FORMAT ('§ INPUT INDEX AND STANDARD DEVIATION") : 
FORMAT (" THE LAST"eI25" INPUTS HAVE BEEN IGNORED") 
FORMAT (' INPUT COLUMN ¢« ROW AND CORRELATION") 
FORMAT (" INPUT MFAN VALUE FOR INDEX ''5]3) 
FORMAT ("§ INPUT STANDAQN DEVIATION FOR INDEX "673) 
FORMAT ('* INPUT ROW AND CORRELATION FOR COLUMN '"'412) 
FORMAT ("© INPUT CORRELATION FOR ROW "¥oIT2e" COLUMN "412) 
END 
SUBROUTINE NAMEFL (IFPROR) 
COMMON/MCS4/CORR (78) sANTM (2912) 9 V (25) sNMEAS oO KNo 

SQRNOT (78) oX (25) »s TUseNTMEAS eo ICOUNT ¢ 

INAME (4) e I TMEASeNIMEASoNTEMP 
COMMON/MISC/IOINs IOOUTs ICARD 


t2 





AYAAaADA 


100 


110 


130 


140 


Cc 
150 


160 
170 
c 

1R0 
190 


AHANAANAAN 


100 


110 


##ePROMPTS USER FOR NAME OF SAMPLE FILE 


DIMENSION IRUFF (4) eXVAL (4) 
DATA JRUFF /4°7]0H / 
IFRROR=0 


#eINITIALIZE CONTENTS OF [8UFF (4) WITH 
#eORL ANK CHARACTERS 


DO 100 J2=164 
IAUFF (J) =10H 
CONTINUE 


CONT INUE 

IF (IOIN.ZEQ.S) WRITE (66180) 

CALL INPT (4e4eNVALSZsI BUFF + XVAL® IEOF ) 
IF (IEOF LT-2.0R-ITEOF.EQ.S5S) GO TO 150 
IF (IOIN.NE.S) GO TO 140 

GO TO (12051300170)%5 IFOF=)] 

CONTINUE 

CALL HFLP (14) 

GO TO 110 

CONTINUE 

WRITE (62190) 

GO TO 110 

CONTINUE 

WPITE (69190) 

IFRROR=)] 


##@LQOAD CONTENTS OF IBUFF INTO INAME ARFAY 


CONTINUE 
NN 160 Jeles 
INAMF (J) =I 3UFF (JU) 


CONTINUE 

RETURN 

FORMAT ('' OPERATOR SAMPLE DESCRIPTION(40 CHAK MAX) ----") 
FORMAT (" INVALID RESPONSE") 

END 


SUBROUTINE OPSAMP (ISAMP) 
#eeGENERATE OPERATOR SAMPLE 
COMMON /MCSM/CORP (78) e ANTM (2912) 0 V (25) 9NMEAS oO KN oe 
1 SQROOT (78) «X (25) sIUeNTMEASe ICOUNT 9 
2 INAME (4) g I TMEASeNIMEASeNTEMP 
COMMON/MISC/IOINe IONUT+ I CARD 
DIMENSION JFRAN(9) 
DATA (ITRANCI) 0 F219) 43333333395 290555 I7T7T7T TTI e2222222 1 
1. 44444444 3 9646666669 9 B8EBBEEBS 2 999999999012 3456789/ 


##01 TNEAR ARRAY CORR(78) IS LOABED INTO LINEAR 
2eeQUMMY ARRAY SQROOT(7&) BEFORE MATRIX TRANSFOR]= 
#@eMATIONS, ALL MATRIX TRANSFORMATIONS ARE 
#eePEFREFORMFI) ON SQROOT(78) TO PRESERVE CORR(78) 
##2AS IMPLEMENTED WITHIN THE EDITING PROCESS. 


IOIN=5 

NO 100 K=1,KN = 
SOROOT (K) =CORR(K) 

CONTINUE 


oo PERFORM MATRIX TRANSFORMATIONS 
CALL SOQPMX (NERROR) 

IF (NERROR.EQ.1) RETURN 

#%% ACCEPT OPERATOR SAMPLE SIZE 
CONTINUF 

WPITE (Ae320) : 
CALL INPT (lelseNVALSZseNSAMPoXVAL » IEOF ) 
GO TO (15024]702¢130214027P90214N)- IFOF } 
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120 
130 
140 


150 


Cc 
C 
c 
160 


170 
180 


190 


200 


210 
220 
250 


240 


250 


260 
C 
C 
c 
c 


270 


260 
G 


CONTINUE 

WRITE (69330) 
GO TO 110 
CONTINUE 

CALL HELP (15) 
GO TO 110 
CONTINUE 
WEITE (69340) 
GO TO 110 
CONT INUF 
IF (NSAMP.G6T.0) GO TO 160 
WEITE (62350) 
GO TO 110 


@@@ASK USFR FOR MACHINE TYPE (32 OR 60 KIT) 


CONTINUE : 
WPITF (69300) 

CALL INPT (lel eNVALSZeMTYPE oXVALe JEOF) 
GO TO (1904160+170.180+290+180)+ IEOF+) 
CONTINUE 

C4LL HFLP (20) 

60 TO 1360 

CONTINUE 

WPITE (60340) 

GO TO 140 

CONT INUF 

IF (MTYPE NE se 32e AND CMTYPENE-460) GO TO 180 
CONTINUE 

WOITE (40310) 

#e% OPTION FOR ACCEPTING A RANDOM SEED. 
CALL INPT (leleNV4LSZeTSoXVAL oILEOF ) 

60 TO (2404210 922052305290+230)+ IEOF?)} 
CONT INUE 

WEITE (62330) 

GO TO 200 

CONTINUE 

CALL HELP (212) 

GO TO 200 

CONT INUE 

WRITE (69340) 

GO TO 200 

CONTINUE 

IF (IS.EFE0-0) GO TO 250 

IF (TSeLTeleO8.1S-GT.9) GO TO 239 
TU=TRAN (IS) 

60 TO 240 

CONT INUF 

1U=67111133 

CONTINUE 

REwIND NTEMP 


feoeWOITTE 40 CHARACTER NAME FIELD INTO 
COCFILE SAMPLE (NTEMP) 


WRITE (NTEMP) INAME eNMEAS 


00 280 J=l+NSAMP = 
#o0 STMULATE A PSEUDO SAMPLE CASE. 

CALL SIMSUS (MTYPE) 

DO 270 K=)l—eNMEAS 

VK) SV (CK) PANTM (20K) eANTM( 1] oK) 

CONTINUE 

WRITE (NTEMP) (VC IO) oe I90=] eNMEAS) 

CONTINUE 


eeeSFT ISAMP TO 1 TO INDICATE OPERATOR SAMPLE HAS BEEN GENERATED 


ISAMP=] 
FND FJIF NTF MP 
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REWIND NTEMP 

REWIND NTMEAS 

WPITE (62360) 
250 CONTINUE 

RE TURN 


300 FORMAT (* ENTER MACHINE WORD SIZE (32 OR 60 BIT) =") 
310 FORMAT (" ENTER RANDOM NO SEED (0 IF NONE OR 19) <=") 
370 FORMAT (¢" INPUT SAMPLE SIZE -!") 
330 FORMAT (* NO INFORMATION ENCOUNTERED") 
340 FORMAT (" INVALID RESPONSE?!) 
350 FORMAT (" SAMPLE SIZE MUST BE GREATEF THAN 0O'*) 
360 FORMAT (" SAMPLE DATA SAVED ON FILE SAMPLE OR TaPe 3.'*) 
END 
SUBROUTINE PRINT 
c 
C *###QOUTPUT OF SIMULATED OPERATORS : 
C 
C #@@NOTE: ALL OPERATOR MEASUREMENTS 
C *##fARE STORED ON FILE SAMPLE. PRESENTLY EACH RECORD 
C #@# CONTAINS NMEA™ ANTHROPOMETRIC MEASUREMENTS 
c 
COMMON/MCS4/CORR (78) sANTM (2912) ov (25) eNMEAS SKN 
l SOROOT (78) eX (25) sIUsNTMEAS es ICOUNT » 
2 INAME (G4) eI TMEASeNIMEASeNTEMP 
COMMON/MISC/IOINe IOOUT + I CARD 
INTEGEP FMT 
DIMENSTON FMT (46) eNUM(12) eTEMP(12) s INDEX (2) 
DATA FMT/4H( 94H 04H (6X) 04H 94H(F6.04H3)) / 
NATA NUM/1H] oe LH2e1H3¢lHGelHSsLH6s1H7 so LHe LHS s2HlOe2H1ls2Hie/s 
ISTOSES0 
IOPLAG=0 


Cc 
Cc #@@ MENU CHOICE OF OUTPUTS 
c 
) 


00 CONTINUE 
WPITE (62630) 
CALL INPT (LeleoNVALSZeTOPT sXVAL » IE OF} 
GO TO (12091109110911053900110)s IEOF?} 
110 CONTINUE 
WRITE (69650) 
GO TO 100 
120 CONTINUE 
IF (IOPT.LT~1-0R.~IOPT.GT.%) GO TO 110 
IF (IOPT-EQ.-1) GO TO 130 
IF (IOPT.E9.2) GO TO 150 
IF (IOPT.-EQ.4) GO TO 390 
GO TO 190 


#@@ PRINTOUT OF MEANS/STANDARD DEVIATIONS 


KOO 


30 CONTINUE 
WPITF (69580) 
DO 140 JROW=] »sNMEAS 
WRITE (620660) IROWs (ANTM(J0sIROW) 9s 021 92) 
140 CONTINUE 
WPITE (62620) 
GO TO 100. 


#¢@ PRINTOUT OF ANTHROPOMETRIC CORRELATIONS 


—“Aagn 


So CONT INUF 
WPITE (69590) 
WRITE (6296550) (NUM(10) -10=] sNMEAS) 
WRITE (69600) 
NO 180 TROW=SteNMFAS _ 
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NO 160 TCOL=]+«NMEAS 
N=NDX (TPROW,ICOL) 
TEMP (ITCOL) =CORRI(N) 
160 CONT INUE 
IF (IROW.GT.1) GO TO 170 
WPITE (69610) (TEMP CICOL) s ICOL=I~OweNMEAS) 
GO TO 1A0 
179 FMT (2) =NUM(TROW=]) 
FMT (4) =NUM(NMEAS¢]<-] ROW) 
WRITE (AeFMT) (TEMP(ICOL) sICOL=I 2OWsNMEAS) 
140 CONTINUE é, 
WPITE (69620) 


GO TO 100 
Cc 
C ##*® PRINT OF ANTHROPOMETRIC MEASUREMENTS 
c 


190 CONT INUF 
WPITF (60430) 
IF (IOFLAG.LT.1) WRITE (69440) 


INOFLAG=] 
Cc 
C #80 TNPUT CARPIT ACT RETURN -- ENTIRE SAMPLE LIST 
C ##e8INPUT SINGLE NUMBER -- SINGLE OPERATOR 
C *8e7NPUT TWO NUMBERS -- SEGMENT OF OPERATOR LIST 
Cc 


CALL INPT (le2@sNVALSZe INDEX »XVAL » IEOF ) 
GO TO (23042709200+210+1009220)s IEOF¢} 
200 CONTINUE 
CALL HFLP (19) 
IOFLAG=0 
GO TO 190 
210 CONTINUE 
WRITE (62420) 
GO TO 190 
220 CONT INUE 
IOLISTZNVALSZ=-2 
WEITE (69450) IOLIST 
230 CONTINUE 
C s#¢eTNDEX VALIDATION 
IF (INDEX(1).GT.0) GO TO 240 
WRITE (69480) 
GO TO 190 
240 CONTINUE 
C 
C *88TO DETERMINE WHETHER FILES WILL BE REWOUND 
Cc 
IF CISTORE.GTeN-ANDeINDEX(1) -GTeISTORE) GO TO 250 
NRSKIP=INDEX (1) =] 
GO TO 260 
250 CONTINUE 
NPSKIP=INDEX (1) -ISTORE-] 
260 CONTINUE 
ISTARTH=INDEX(1) 
C **0ITF TWO INPUTS ENCOUNTERED = OUTPUT LIMITS 
C *#88TF ONE INPUTS ENCOUNTERED = OUTPUT SINGLE OPERATOX 
IF (NV4LSZ.EQ.1) GO TO 280 
IF (NVALSZ.GE.2) GO TO 290 
270 CONTINUE 
ISTORE=0 
ISTART=] 
IEND=999 
NRSKIP=0 
GO TO 320 
280 CONTINUE 
IEND=ISTART 
an TO 30 
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290 CONTINUE 
IF (INNEX(1) LES INDEX(2)) GO TO 300 
WPITE (62490) 
GO TO 190 
300 CONTINUE 
IEND=INDEX (2) 
310 CONT INUF 
IF (INNEX(1) -GTeISTORE. AND. ISTOREeNE-0) GO TO 330 
320 CONTINUF 
C #8#8REWIND FILE. 
REWIND NTEMP 
READ (NTEMP) INAME sNMEAS 
IF (EOF (NTEMP) .EQ.1) GO TO 410 
330 CONT INUE 
IF (NPSKIP.LT.1) GO TO 350 
C #8°9WIND FILE FORWARD 
00 340 Y=l]1eNRSKIP 
READ (NTEMP) (V(K) eK=] «NMEAS) 
IF (FOF (NTEMP) .EQ.1) GO TO 400 
340 CONTINUE 
350 CONTINUF 
C #e#8f20UTPUT HEADER 
WRITE (62530) (INAME (IO) eI0=]1 64) 
IF (ISTART.EQ.LIEND) WRITE (62540) IEND 
IF (ISTART.NE.TEND.-ANDJLIEND.NE 999) WRITE (64550) ISTARTsIEND 
WPITF (62500) . 
WRITE. (629660) (NUM(IO) «- [O=]) ¢eNMEAS) 
WRITE (62600) 


c 
C #8#80UTPUT NATA FROM FILE 
S 

ISKIP=9 

NFOUND=0 


99 360 Y=ISTART+«IEND 
C #88#SKIP LINE FEATURE 
IF (ISKIP-£Q.1) WRITE (6460) 
ISKIP=0 
IF (ISTART.NEJTEND.AND.I NE SIENDe ANDY (IT /5)%S5.EQ21) ISKIP=) 
READ (NTEM®) (V(K) «K=] eNMEAS) 
IF (FOF (NTEMP) .FQ.1) GO TO 370 
NFOUND=NF OUND +] 
ISTORE=I 
WOTTE (69510) (V(K) oK=] eNMEAS) 
350 CONTINUE 
GO TO 3A0 
370 CONTINUE 
IF (TEND eNE.}S99DCANDeISTART.LNELIEND) WRITE (69470) IEND 
ISTORE=I 
380 CONTINUE 
WRITE (62520) 
GO TO 190 
390 CONTINUE 
REWIND NTEMP 
RF TURN 
400 CONTINUE 
IOLIST=I~2e¢I STORE 
WPITE (60560) 
GO TO 190 
410 CONT INUE 
WRITE (69570) 
GO TO 100 


420 FORMAT ("* INVALID RESPONSE") 


430 FORMAT (" OQUTPUT MODE: <<") 
440 FORMAT (" FOR ALL OPFRATOR NATA (ENTER CARRIAGE RETURN) "o/9 


ad 





450 
460 
470 
480 
490 
500 
510 


530 
540 
550 
560 
570 


540 
590 


600 
610 
670 
430 
640 
650 
660 


AAAAN 


100 


Aa- 


AAO 


]* 


FORMAT 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
END 


SUBROUTINE RAND2 


FOR INDIVIDUAL OPERATOR DATA 
2" FOR A SEQUENCE OF OPERATORS 


(ENTER OPERATOR NUMBER) M9/ 
(ENTER RANGE OF OPERATOR NUMBERS) "/ 


* " END (ENTER $)") 
FORMAT (* THE LAST'tsI46% INPUTS RAVE BEEN IGNORED") 
FORMAT (1K ) 
FORMAT (" ERROR IN UPPER INDEX |= 414) 
FORMAT ('§ INDEX NUMBER CANNOT BF LESS THAN OR EQUAL TO 0") 
FORMAT ('* FIRST INOEX MUST BE LESS THAN SECOND INUVEX") 
FORMAT (27Xe"OPERATOR SAMPLE''o/s71(1LH=) 9) 
FORMAT (1lIFG.29F 5.2) 
FORMAT (71(1H=)) 
FOPMAT (1% eGA10010XeA10) 
FORMAT (l10Xe*OPERATOR NO» 4) 
FORMAT (1OXe"OPERATOR NOS eIGq" w=, 14) 
FORMAT ("§ INDEX NUMBER EXCEEDS RECORDS IN OPERATOK FILE") 
FORMAT ('* FILE SAMPLE IS NOT ATTACHED TO THE PROGRKAM!!/ 
" OPERATOR OUTPUT PEQUEST ARORTED --="9//) 
FORMAT (/e21])Xe"MFEANS AND STANDARD DEVIATIONS'*e/91Xo70(1H=) 9/92X 9M 


LE ASUREMENT'* 6X op "MEANS" 6X e"STA DEV oo 1X0 71 (LH=) ) 


(7e7Xs"CORRELATION MATRIX FOR SO RO Es he MEASUREMENTS" 97 


LelXe70(1H=)) 


(1Xe70(1H-) ) 

(l12F6.3) 

(1X70 1H=) ) 

(" PRINT MONE (1-MEAN/STD32-CORR ato OPER 
(I1995Xe2F12.2)) 

(" INVALIN RESPONSE?) 

(2X e911 (A2)94X) sA2) 


MEAS34=-END) <'') 


(RNORM oMTYPE «QUAN, I J) 


@@eGFNERATES A NORMAL DISTRIBUTION wITH MEAN OF ZERO AND 
STANDARD DEVIATION EQUAL TO ONEcece. 
®eeMETHODOLOGY --= MARSAGLIA = BRAY 


COMMON/MCS4/CORR (78) pANTH (2912) 0 V¥ (25) oNMEAS OKNe 

l SQROOT (78) eX (25) eJUsNTMEAS oe ICOUNT » 
2 INAME (4) sITME AS »NIME AS eNTEMP 

DIMENSTON QUAN (2) 

COMMON/MISC/IOINs IOOUT+ICARD 

TJ=I Je] 

IF (LJ.FEQ@e2) GO TO 110 

CONT INUF 

X1=2.0@UNFORM(MTYPE) -1.0 

YIS2.N*UNF ORM(MTYPE) -1.0 

S=X1*#X1]+Y127Y) 

IF (SeGFe1.0) GO TO 100 

S=SQPRT (~2.07AL0G(S)/S) 

QUAN (1) =X145S 

QUAN(2)=Y14S 

CONT INUE 

RNORMSQUAN (TJ) 

IF (IJeFQe2) JU=0 

RETURN 

END 

SUBROUTINE REASON (CICALLeICOLsI]ROWsTSTOREL] + TSTORECs IFATAL) 
@e%e CHECKS FOR REASONABLENESS OF INPUT VALUES 
COMMON/MCSM/CORR (79) eANTM( 2912) 9 V(25) sNMEASSKNG 
] SQROOT (7A) oX (25) eIUeNTMEASe ICOUNT 9 
2 INAME (4) s I TME AS oNIMEAS oNTEMP 
COMMON/MISC/IOINs IOOUTs ICARD 


1 MEANS AND STANDARD DEVIATIONS 
CORREI ATTON MATPRLX. VERPTFICATION 


ICALL = 
TCALI = ? 
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c 


CY CY Or -— OOM KH OOANMNAN OaAaN 


mHOOoM 


20 


130 


140 


150 


160 


ICALL = 3 INDEX VERIFICATION 


IFATAL=0 

IAFLAG=0 

IF (ICALtetLT.0) IAFLAG=] 
ICALL=JABS(ICALL) 

IF (ICALLeGT.1) GO TO 110 


2 MEAN VALUES MUST BE WITHIN 4 THREE STANDARD VEVIATION 
ee NETGHSORHOOD OF THE ACCEPTED MEAN VALUES 


IWRITE=] 

IF (TROW.LT.1eOReIROW.GTeNMEAS) GO TO 120 
IF (ICOL-EO.2) GO TO 100 

RETURN 


#@% STANDARD DEVIATIONS MAY NOT BE NEGATIVE 
#9 ANT THE MEANS MUST BE GREATER TRAN 
#00 THE STANDARD DEVIATION. 


CONTINUE 

IWRITE=2 

IF (TSTORE2?.LT.0-) GO TO 120 
IWRITE=3 

IF (TSTORE].LE.TSTORE2) GO TO 120 
RE TURN 


#@@ CHECKS VALIDITY OF INDEX NUMGERS 


CONTINUE 
IWRITE=] 
IF (IROWeLE.0-OR-IROW.GTeNMEAS) GO TO 130 
IWRITE=? 
IF (ICOLeLE/0-0PeICOL«GTeNMEAS) GO TO 130 


#eeTS THE NUMBER WITHIN THE UPPER TRIANGULAR REGION? 


IWRITE=3 
IF (IPN04.GT.ICOL) GO TO 130 
IF (ICALL-E&O.3) RETURN 


#0 ALL DIAGONAL ELEMENTS MUST EQUAL 1.000 
#o0@ MATRIX ELEMENTS MAY NOT BE LESS THAN (-1) OR GREATER 
oo THAN ONE. 


IWRITEs4G 

IF (ICOLeEQeIROW. AND. TSTORE].NE.1-000) GO TO 130 
IWRI TESS 

IF (TSTORE]LTe~-1.000) GO TO 130 

IWRITE=6 

IF (TSTORE]-GT.1.000) GO TO 130 

RF TURN 


#f#eRRANCH TO APPROPRIATE ERROR MESSAGE 


CONTINUF 

IF CIWRITE.NE) eAND-IOINONE.S) WRITE (62280) [ROW 
GO TO (19041404150) s JIWRITE 

CONTINUE 

GO TO (1904200921051600170+180)s IWRITE 
CONTINUE 

WPITE (42230) 

GO TO 220 

CONTINUE 

WRITE (69240) 

GO TO 220 

CONTINUE 

WEITTF (He2SN) TROW-TCHI 
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170 


130 


190 


200 


210 


220 


230 
240 
250 


250 
270 
280 
290 


300 
319 


aAaANDN 


GO TO 220 

CONT INUE 

WRITE (69260) IROWeICOL 
GO TO 220 

CONTINUE 

WPITE (69270) IROWsICOL 
GO TO 220 i ' 
CONTINUE 

WRITE (6290) JROWsNMEAS 
GO TO 220 

CONT INUF 

WEITE (62300) ICOL»NMEAS 
GO TO 270 

CONTINUE 

WEITE (69310) TIROW,sICOL 
GO TO 220 

CONTINUE 

IFATAL=] 

RETURN 


FORMAT (" STANDARD DEVIATIONS LESS THAN ZERO!) 

FORMAT ("§ STANDARD DEVIATION GREATER THAN MEAN VALUE."') 
FORMAT (*§ DIAGONAL ELEMENT IN ROW "oI25" COLUMN "sI2¢ 

+ "' [TS NOT FQUAL TO 1.0") ; 

FORMAT (" ENTRY IN COLUMN "oI2e" ROW "sles 


Seems cos THaN —),0'*) 
FORMAT (" ENTRY IN COLUMN "eI2s" ROW "sI2s 
« " TS GREATER THAN #1,0') 


FORMAT (110) 
FORMAT ("* 20W "sI26" IS LESS THAN ZERO OR GREATER THAN "912) 


FORMAT (" COLUMN "“,T2e" IS LESS THAN ZEPO OR GREATER THAN "97 2) 
FORMAT (" ROW "eT2o' COLUMN "sI2e" IS OUTSINE OF UPPER", 

+ ™ TRIANGULAR REGION OF CORRELATION MATRIX." ) 

END 


SUBROUTINE SAVE 
#@0 THIS ROUTINE SAVES INPUT DATA ON PERMANENT FILE TAPE? 


COMMON/MCSM/CORP (78) «ANTM (20¢12) 9 V (25) sNMEAS oKNo 
1] SQROOT (78) 9X (25) sIUsNTMEASs ICOUNTs 
2 INAME (4) sITMEASsNIMEAS so NTEMP 
COMMON/MISC/IOINs IOOUTs I CARD 
REWIND NTMEAS 
@OOSAVE NO# OF MEAS.e FLEMENTS IN CORR MATRIXe SEX FLAG 
WRITE (NTMEAS) NMEASeKNo INAME 
#80 SAVE MEAN/SSTD. VECTOR, 
WRITE (NTMEAS) ((ANTM(I0¢J0) sJO=1 92) sJO=1 sNMEAS) 


C ##% SAVE CORRELATION MATRIX 


mm C) 
QO 
So 


ANAINANANANRAN 


WRITE (NTMEAS) (CORK(IO) » IO0=1 6KN) 
REwIND NTMEAS 

WRITE (69100) 

RETURN 


FORMAT ( MEASUREMENT DATA SAVED ON FILE UMEAS OR TAPE2") 
END 
SUBROUTINE SIMSUB (MTYPE) 


oe 


#@3 FOR A GIVEN SQUAR ROOT MATRIX R OF ORDER M IN UPPER 
#@0 TRIANGULAR FORM OF A CORRELATION OR A COVARIANCE 

#o% MATRIX OF M COLUMNSs THIS ROUTINE COMPUTES A VECTOR 
#o0 6V OF PSEUDO SCORES OF LENGTH M SUCH THAT IT WILL 

#00 BE A QUAST=RANDOM S4MPLE FROM A NORMAL POPULATION 
#20 WITH A MEAN VECTOR OF ZEROS AND THE SAME CORRELATION 
#@% OR COVARIANCE MATRIXs ITeEe FROM AN N(OoR) MAIRIX.S 


COMMON/MCS4/CORR I 7R) oe ANTM (26127) 0 V (25) sNMFASeKNe . 
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100 


110 


AANAAADND 


AAAN 


100 
c 


120 
130 
e 


140 


} SQROOT (78) eX (25) sIUseNTMEAS » ICOUNT >» 
2 INAME (G6) s I TMEAS eNIMEAS eNTEMP 
COMMON/MISC/IOIN»s IOOUT s I CARD 
OIMENSION Y(45) 

IJ=0 

DO 100 I=leNMEAS 

V(I?=0.0 

CALL RAND2 (RNORMeMTYPE ,QUANes IJ) 
Y (I) =RNORM 

D0 110 [=leNMEAS 

90 110 J=lel 

K=NDX (Jel) 

VIII HV(I) ©¥ (I) @SOROOT (K) 

RETURN 

END 

SUBROUTINE SQRMX (NERROR) 


eee FOR A GIVEN SYMMETRIC MATRIX OF ORDER N STORED IN ITS 
#ee UPPER TRIANGULAR FORM. COLUMNWISE IN THE VECTOR R OF 
@@0 LENGTH N®(NO1)/25 THIS ROUTINE TRANSFORMS IT INTO ITS 
eee MSQUARE ROOT" IN THE SAME FORM IN THE veCTOR SQROOT(78) 


COMMON/MCSM/CORR (78) eANTM (2012) 0 V(25) eNMEASOKNe 
1 SQROOT (78) «X (25) sIUeNTMEAS es I COUNT 
Z INAME (4) es I TMEASeNIMEAS eNTEMP 
COMMON/MISC/IOINs IOOUTs I CARD 
DIMENSION ZX(7) 
INTEGER ZX 
@eONFRRORP FLAG IS USED TO SIGNIFY A CONDITION 
@OOWHICH WILL PRODUCE A NEGATIVE SQRUOTI(N) VALUE 
feo JUST BEFORE IMPLEMENTATION OF THE SQUARE ROOT 
SOOFUNCTION. 
NF RROR=0 
DO 120 [=leNMEAS 
ZX (G4) =7T 
N=NDX (I eI) 
SeOTEST FOR NEGATIVE PANDICAL ARGUMENT 
IF (SOQROOT(N).LT.9) GO TO 130 
SQ=SORT (SQROOT(N) ) 
DO 100 J=IeNMEAS 
N=NDX(T+J) 
SQROOT (N) =SQROOT (N) /SQ 
CONTINUE 
SOONORMAL FXIT FROM SUBROUTINE 
IF (ITeFQOeNMEAS) NERROR=O0 
IF (IL eEQSNWEAS) RETURN 
ITl=I +} 
ZX(5)=T2 
NO 110 J2=I1leNMEAS 
00 110 J2=l[2@eNMEAS 
N=NDX (J 2eJ2) 
L=NDX (I eI2) 
K=NDX (T »J2) 
paceorin jesaroanintesocooriesoRooTiK) 
ZX(6)= 12 
ZX(7)=J2 
CONTINUE 
ZX(1)=N 
ZX(2) EL 
2X (3) =K 
CONTINUE 
CONTINUE 
SOOARNORMAL EXIT FROM SUBROUTINE 
NF RROR=)] 
WRITE (60190) 
CONTINUE 
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150 


160 


170 


180 
190 


200 


210 
220 


100 


WRITE (49190) J 
CALL INPT (Goel eNVALSZeT ANS eXVAL eo IEOF ) 
GO TO (1704160915041 70+1605170). LEOFe] 
CONTINUE 

WPITE (69200) 

WRITE (69210) (ZX( 1) 9IT=107) sSQROVUTI(N) 


CONTINUE 

RETURN 

CONT 1NUE 

WRITE (69220) 

GO TO 140 

FORMAT ('§ ANTHROPOMETRIC CORRELATION VALUES ARE INVALID") 

FORMAT (" #=£ ENTER CARRIAGE RETURN TG CONTINUE PRUCESSING"9/ 
Wome ENTER ? TO DISPLAY SQRMA VARIABLES") 

FORMAT ("§ AN INDEFINITE QUANTITY HAS BEEN DETECTEU BY SQRMX"/ 
" VARTABLE LIST WITHIN SORMX:"/" N Lit. 

on K I Il le J2 SQOROOT" esol OX e44(1]1H=) ) 


FORMAT (BX e415 ol Xe STS oF Ge 3 eS 910K 944( lH) s//) 
FORMAT (/" INVALIN RESPONSE FOR" 7) 
ENO 

FUNCTION NOX(T,J) 

IF (I1-e6TeJ) GO TO 100 
NOA=T + (J®# =) /2 

RF TURN 

NDXsJe(T#Il-1) 72 

RF TURN 

END 

FUNCTION UNFORM(MTYPE) 


SeeUNTFORM RANDOM NUMBER GENERATOR ON INTERVAL (0091) 


COMMON /MCSM/ CORR (78) eANTM (2912) 9V(25) sNMEASoKNo 


1SQRGOT (78) oX (25) sIUeNTMEASe ICOUNT 9 
CINAMF (G4) el TMEASeNIMEAS sNTEMP 


COMMON/MISC/IOINs LOOUT+ICARD 

60 BIT UNIFORM RANDOM NUMBER GENERATOR 
DATA I1X/16777213/9C/281474976710655./ 
IF (MTYPE.EQ.32) GO TO 100 

[U=IU#]IX 

UNF ORM=TU/C 
RF TUPN 

32 BIT UNIFORM RANDOM NUMBER GENERATOR 
[U=IU#65539 

IF (IU.GEe0) GO TO 110 
TU=1U¢2147483647¢1 

UNF ORM=TU 

UNFORMZUNF ORM®. 465661] 3F -09 

RETURN ee oo 

END 

SUBROUTINE HELP (T) 

RE TURN 

ENO# 
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oon ouvuckwhn = 


100 


110 


120 


130 


140 


150 


160 


170 


180 


190 


200 


APPENDIX D 


Program Listing for CAPE Model 
Modified for Use with FORTRAN IV 


COMMON/MCSM/CGRR (78) , ANTM (2,12), V (25) , NMEAS, KN, 


1 SQROGT (78) ,X (25) , 1U, NTMEAS, ICOUNT, 
2 INAME (10), 1TMEAS, NIMEAS, MAXSAM,NTEMP 


COMMON/MISC/IOIN, 100UT, ICARD 

DATA NUCREW/8/ 

DATA MAXFUN/Y/ 

DATA MAXOPT,MAXANT, MAXACT/3, ?,3/ 
100UT=0 

ISAMP=0 

MMEAS = 12 

CONTINUE 

IGIN=5 

WRITE (6,290) 

CALL INPT (1,1,NVALSZ,NMEAS, XVAL, JEOF) 
IEOFPI=l1EOF+1 

GO TO (130,100,110,120,280,120), IEOFP! 
CONTINUE 

CALL HELP (1) 

GO TO 100 

CONTINUE 

WRITE (6,340) 

GO TO 100 

CONTINUE 

IF (NMEAS.LT.1.QR.NMEAS.GT.MMEAS) GO TO 120 
KN=NMEAS™ (NMEAS+1) /2 

1INPT=0 

CONTINUE 

IGIN=S5 

CONTINUE 

WRITE (6,300) 

CALL INPT (1,1,NVALSZ,10PT, XVAL, 1EOF) 
IEGFPI=1E0F+1 

GO TG (180,140,160,170,280,170), 1EOQFP! 
CONTINUE 

CALL HELP (2) 

GO Td 150 

CONTINUE 

WRITE (6, 340) 

GO TO 150 

CONTINUE 

IF (IOPT.LT.1.0R.10PT.GT.MAXANT)GO TO 170 
GO TO (190,200,220, 240, 260, 280, 225) , 1OPT 
CONTINUE 

CALL ACCPT (ISAMP) 

TINPT=1 

GO TC 140 

CONTINUE 

IF C(LINPT.EQ.1) GO TO 210 

WRITE (6,350) 

GO.TO 140 
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a 
32 
93 
S4 
95 


97 
58 
99 
60 
61 
62 
63 
64 
65 


67 
68 
69 
70 
71 
72 
73 
74 
75 


77? 
78 
79 
80 
81 
B2 
83 
BY 
B85 
86 
87? 
88 
89 
90 
91 
92 
33 
94 


210 


220 


e2e 


225 


230 


250 


260 


270 


280 


290 


300 


310 
320 


330 
340 
350 


360 


CONTINUE 

CALL EDIT 

GO TO 140 

CONTINUE 

Pel oOANee be bieGo 10 230 
WRITE (6,320) 

GO TO 140 

CONTINUE 

IF (ISAMP.NE.1) GO TO 222 
CALL PUNRT 

GO T0 140 

CONTINUE 

CALL SAVE 

GO 10 140 

CONTINUE 

IF CIINPT.EQ.1) GO TO 250 
WRITE (6,360) 

GO TO 140 

CONTINUE 

CALL PRINT 

GO 10 1480 

CONTINUE 

IF (ISAMP.NE.1) GO TO 270 
CALL OPSAMP (ISAMP) 

GO TO 140 

CONTINUE 

WRITE (6,330) 

GO TO 140 

CONTINUE 

WRITE (6,319) 

STOP 

FORMAT (* ENTER NUMBER GF MEASUREMENTS 


FORMAT (° OPTIONS (1-INPUT;2-EOIT; 3-SAVE;4Y-PRINT;5-GENERATE; * 


Peo seNOs 7-PUNEH) °) 
FORMAT (5X, ‘END OF OPERATOR SAMPLE‘) 


FORMAT (1X, *«xMCSM1YQ UNABLE TO SAVE MEASUREMENT GATA - FATAL 


1 ERRORS") 


FORMAT (10X, ‘xx»MCSMI165 UNABLE TO GENERATE SAMPLE - FATAL ERRORS *) 


FORMAT (/, *xx»MCSMOCU INVALIO RESPONSE *) 


FORMAT (10X, *xx*MCSM130 ANTHROPOMETRIC DATA MUST BE INPUT BEFORE™, 


ieox, (HE EQIT FUNCTION. *) 


FORMAT (10X, *xxxMCSM1SO ANTHROPOMETRIC GATA MUST BE INPUT BEFORE", 


1/,13X, ‘THE PRINT FUNCTION. ‘%) 
END 
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oan onuf wn = 


100 


110 


109 


113 


120 


130 


140 


150 


160 


SUBROUTINE INPT (ITYPE, IVALSZ,NVALSZ,IVAL,XVAL, JEGF) 
COMHON/MISC/IGIN, 166UT, ICARD 
DIMENSION ICHAR (80) , IVAL (IVALSZ) ,XVAL (IVALSZ) 


DIMENSION NVAL (20) 


LOGICAL » 1 NCHAR (80) 
EQUIVALENCE (NVAL (1) ,NCHAR (1) ) 


DATR IBLK, IQUES, IDOL, 1QUG,I1POST, IPO, ICOM, IMIN, IPLUS/IH ,1H?, 


L1HS,1H*, 1H: 1H.,1H,,1H-,1Ht/ 


DATA IN8LK/4H / 

ISTBP=0 

IF (IVALSZ.GE.0) GG Td 100 
ISTOP=1 

IVALSZ=IABS (IVALSZ) 

CONTINUE 

1QU8TE=0 

NVALSZ=0 

IEGF=0 

IF (100UT.EQ.1) GO TO 120 

D6 110 I=1,80 

ICHAR (1) =I BLK 

CONTINUE 

READ (IGIN, 330,END=111) ICHAR, NCHAR 
DG 109 J=1,80 

IF (JCHAR (J) .EQ. IBLK} GO TO 109 
G6 16 130 

CONT INUE 

G6 16 113 

CONT INUE 

REWINO IOIN 

JEOF = 1 

RE TURN 

CONT INUE 

1121 CARO 

GO 710 140 

CONTINUE 

111 

CONTINUE 

12=80 

CONTINUE 

WRITE (6,501) IOEC,ISTOP, IVALSZ,NVALSZ, I1,132, IEOF 
IF (ISTOP.EQ.1.AND.IVALSZ.EQ.NVALSZ) RETURN 
IF (11.LE.12) GO TO 160 

IF (NVALSZ.LE.IVALSZ) RETURN 
IEOF=5 

RETURN 

CONT INUE 

DG 170 I=!I1,12 
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170 


180 


190 


200 


210 


220 


230 


240 


250 


IF ¢€ICHAR(]) .EQ.IBLK ) GO TO 170 
IFIRST=] 

GO TO 1860 

CONTINUE 

IF (NVALSZ.EQ.0) IEOF=1 

IF (NVALSZ.GT.IVALSZ) JEOF=5 

RETURN 

CONTINUE 

IF €ICHAR (IFIRST) .EQ.IQUES) JEOQF=2 
IF (ICHARCIFIRST) .EQ.IDGL) JEOF=4 
IF (IEQF.NE.0) RETURN 

IF (ITYPE.EQ.4) GO TO 300 

IF (ITYPE.NE.3) GO TO 190 

IF ¢ICHAR CIFIRST) .NE.1QU0.AND.ICHAR (IFIRST) .NE.IPOST) GO TO 190 
IFIRST=IFIRST+1 

JQUOTE=1 

CONTINUE 

DG 220 J=IFIRST,1e 

IF (1QUOTE.EQ.0) GO TO 200 

IF ¢€ICHAR (1) .€0.]QU0.0R. ICHAR (1) .EQ.IP0ST)GO TO 210 
GO TG 220 

CONTINUE 

IF (ICHAR (1) .NE.IBLK.AND.ICHAR (1) .NE.1CGM) GO TO 220 
CONTINUE 

ILAST=1-1 

GO TO 230 

CONTINUE 

ILAST=12 

CONTINUE 

NVALSZ=NVALSZ+1 

TizILAST+2 

IF (NVALSZ.GT.IVALSZ) GO TG 150 

IF (ITYPE .€Q.3) GO TO 310 

ISIGN=1 

IF (ICHARCIFIRST) .NE.IMIN. AND. ICHAR (IFIRST) .NE.IPLUS) GO TO 240 
IF (ICHARCIFIRST) .EQ.IMIN) ISIGN=-1 
IFIRST=IFIRST+1 

CONTINUE 

INUM=0 

JOEC=0 

IOECF=0 

06 260 J=IFIRST, ILAST 

IF €ICHAR (1) .NE.IP0) GO TO 250 
JOECF=IDECF +1 

1OEC=1 

GO TG 260 

CONTINUE 

CALL ICONV (CICHAR (1), JADD, IEOF) 

IF (IEOF.EQ.3) RETURN 
]NUM=INUM10+] ADO 
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101 
102 
103 
104 
105 
106 
10? 
108 
109 
110 
111 
112 
113 
114 
115 
116 
11? 
118 
119 
120 
121 
12e 
123 
124 
125 
126 
127 
128 
129 
130 
$31 
132 
133 
134 
135 
136 
137 
138 
139 
140 
iy} 
142 
143 
144 
YS 
146 
147 
148 
149 


$01 


260 


270 


c 


280 


cc 


290 


300 


909 
910 
310 


320 


904 


311 


313 


330 
340 


FORMAT (/ * TEST *, 10110) 

IF (1OEC.NE.0) IOEC=IDEC*1 

CONTINUE 

IF (IOECF.LT.2) G@ TO 280 

CONTINUE 

WRITE (6,501) IQEC, IOECF 

IEOF=3 

RETURN 

CONTINUE 

INUM=INUM*ISIGN 

XNUM=FLQAT (1 NUM) 

WRITE (6,501) IQEC, IDECF, INUM, IAOO, ISIGN, ITYPE 

IF ((ITYPE.EQ.i.ANO. 1OEC.LE.1).0R. (ITYPE.EQ.2)) GO Ta 290 

1EQF=3 

RETURN 

CONTINUE 

IF (IOEC.NE.0) XNUM=XNUM/10. Ox» (JOEC-1) 

IF (JTYPE.EQ.1) IVAL (NVALSZ) =1F 1X (XNUM) 

IF (ITYPE.EQ.2) XVAL (NVALSZ) =XNUM 

ICARO=ILAST+! 

WRITE (6,501) IOEC,1OECF,INUM,IADO, ISIGN,ITYPE,NVALSZ, 
IVAL (NVALSZ), ICARO 

GO To 150 

CONTINUE 

ILAST=IVALSZ™4 

FORMAT { 2X, 80Al) 

FORMAT (2x, 20AU) 

CONTINUE 

DO 320 J=1, IVALSZ 

IVAL (J) =1UBLK 

CONTINUE 

L=ILAST-IFIRST+1 

NVALSZ=2 (L-1) /1041 

IF (L.GT.4xIVALSZ) GO T@ 270 

ENCOOE (L,340,1VAL (1) ) (ICHAR (I), IFIRST, ILAST) 

j= 0 

FORMAT(* ENCOOE *, 10110) 

00 311 I = IFIRST, ILAST 

jase} 

NCHAR (J) = NCHAR (1) 

CONTINUE 

DO 313 1 = 1,1VALSZ 

IVAL (1) = NVAL(D 

CONTINUE 

ICARO=11 

RETURN 

FORMAT (80A1, T1, 80Al) 

FORMAT (80A1) 

END 
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SUBROUTINE ICONV (JCHAR,11, 1EOF) 
DIMENSION JNUM (10) 

DATA JNUM/1HO,1H1,1H2,1H3, 1HY,1H5,1H6,1H7, 1H8,1H9/ 
FORMAT (//° RAT. ICONV °, A2,2110 /) 
00 100 J=1,10 

IF (JCHAR.NE.JNUM(JJ} GO TO 100 
Ti=J-1 

RE TURN 

CONTINUE 

IEOF = 3 

WRITE (6,11) JCHAR, 11, IEOF 

RETURN 

END 
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100 


110 


120 


130 


140 


150 


160 


170 


180 


190 
200 


90! 


SUBROUTINE ACCPT (ISAMP) 

COMMON/MCSM/CORR (78) , ANTM (2,12), V (25) , NMEAS, KN, 
SQROOT (78) ,X (25), IU, NTMEAS, ICQUNT, 
INAME (10), ITMEAS, NIMEAS, MAXSAM, NTEMP 

COMMON/MISC/IGIN, 10GUT, ICARO 

DIMENSION BUFF (3), TEMP (12) 

DATA BUFF/3%0./ 

ISAMP=0 

ICOUNT=0 

INLINE=0 

1OIN=S 

CONTINUE 

WRITE (6,610) 

CALL INPT (1,1,NOINPT, IMQDE, XINPT, JEQF) 

IEQFPI=IEQF+1 

GO TO (130,110,120,110,490,110), JEOFP! 

WRITE (8,650) 

GO TO 100 

CALL HELP (4) 

GQ TO 100 

IF (IMODE.EQ.1) GO TO 200 

IF (IMOQOE.NE.2) GO TO 110 

CONTINUE 

WRITE (6,620) 

CALL INPT (1,1,NOINPT, JFFORM, XINPT, JEGF) 

IEQFP1<IEOF+1 

GO TO (170,150,160,150.490,150), IEOFP! 

WRITE (6,650) 

GO TO 140 

CALL HELP (5) 

GO TO 140 

IF (IFFORM.EQ.1) GO TO 190 

IF (IFFORM.NE.2) GO TO 150 

REWIND NTMEAS 

READ (NTMEAS,END=180) NMEAS,KN, INAME 

READ (NTMEAS) (ANTM (IG, JO), 10=1,2), JO=1, NMEAS) 

READ (NTMEAS) (CORA(10),10=1,KN) 

REWIND NTMEAS 

WRITE (6,690) 

1SAMP=1 

RETURN 

CONTINUE 

WRITE (6,670) 

GO TO 490 

1OIN=1 

CONTINUE 

CALL NAMEFL (IERROR) 

ICOUNT=ICQUNT+I ERROR 

WRITE(6,901) IMODE, IERROR, ICOUNT 

FORMAT (/ * ACCPT *, 10110) 
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210 


220 
230 


250 


260 


270 


280 


290 


300 


INLINE=INLINE*+1 

IF (IMODE.EQ.1) WRITE (6,500) 

IF (IMGOE.NE.1) WRITE (6,510) INAME 

MAXINP=2 

WRITE (6,90!) IMODE, IERROR, ICOUNT, INLINE, MAXINP, IEOF 
DO 320 IROW=1,NMEAS 

CONTINUE 

IF (IMODE.£Q.1) WRITE (6,530) IROW 

NAMT=0 

CONTINUE 

1TOTIN=MAXINP-NAMT 

ICOL=NAMT +1 

CALL INPT (2,1 TOTIN, NOINPT,IVAL, ANTM (ICOL, IRGW) , IEOF) 
IF (IlEOF.EQ.S5.OR.1EOF.EQ.0) NAMT=NAMT+NOINPT 
WRITE (6,901) IMODE, ITOTIN, MAXINP,NAMT, IEOF, ICOL, ICOL 
IF (IMODE.EQ.1) GO TO 240 

INLINE=INLINE*] 

IF (NAMT.EQ.MAXINP) GO TO 290 

IF (IEOF.EQ.1) GO TO 480 

WRITE (6,540) IROW 

IF (1E0F.E0.5) GO TO 270 

WRITE (6,700) INLINE 

ICOUNT=ICOUNT+1 

GO 170 320 

CONTINUE 

IEQFPI=IEOF+! 

GO TO (280,210,250,260,490,270), IEOFP! 
CONTINUE 

CALL HELP (6) 

GO 10 210 

CONTINUE 

WRITE (6,650) 

GO TO 210 

CONTINUE 

IOLIST=NOINPT-ITOTIN 

WRITE (6,600) IOLIST 

IF (]MODE.NE.1) ICOUNT=ICOUNT*+! 

GO 70 290 

CONTINUE 

IF (NAMT.GE.2) GO TG 290 

WRITE (6,550) 

GO TO 230 

CONTINUE 

IF (NAMT.GT.2) NAMT=e2 

IRERR=0 

060 300 ICOL=1,NAMT 

CALL REASON (1, 1COL,1ROW,ANTM (1, IROW) ,ANTM (2, IROW), IFATAL) 
IRERR=IRERR+ I FATAL 

CONTINUE 

IF (IRERR.EQ.0) GO TO 320 
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101 
10e 
103 
104 
105 
106 
107 
108 
109 
110 
111 
lle 
113 
114 
115 
116 
ily 
118 
119 
120 
121 
l2e 
123 
124 
125 
126 
127 
128 
129 
130 
131 
132 
133 
134 
35 
136 
137 
138 
139 
140 
141 
i4e 
143 
uy 
145 
146 
147 
148 
149 
150 


310 
320 


330 


340 


350 


360 


370 


380 


3390 


400 


IF (IMODE.NE.1) GO TO 310 
WRITE (6,630) IROW 


GO TO 220 

ICOUNT=ICOUNT+IRERR 

CONTINUE 

IF (IMOOE.EQ.1) WRITE (6,560) 

IF (IMOOE.NE.1) WRITE (6,520) 

DO 4Y60 IROW=1,NMEAS 

CONTINUE 

NAMT=0 

MAX INP=NMEAS+1-IRGW 

CONTINUE 

IF (IMOOE.E0.1.ANO. IROW.EQ.NMEAS) WRITE (6,580) IROW,MAXINP,NMEAS 
IF (IMODE.EQ.1.AND.1TROW.NE.NMEAS) WRITE (6,570) IROW,MAXINP, IROW, 
INMEAS 

CONTINUE 

TTOTIN=MAXINP-NAMT 

ISTART=1ROW+NAMT 

CALL INPT (2, I1TOTIN,NOINPT, IVAL, TEMP (ISTART) , IEGF) 
IF (IEGF.EQ.0.0R.1EOF.EQ.5) NAMT=NAMT*NOINPT 
IF (IMODE.EQ.1) GO TO 360 
INLINE=INLINE*1 

IF (NAMT.EQ.MAXINP) GO TO 410 

IF (1E€O0F.EQ.0) GO TO 350 

IF (1EGF.EQ.1) WRITE (6,710) 

IF (IEGF.EQ.1) GO TO 470 

IF (I€0F.EQ.5) GO TO 390 
ICOUNT=ICOUNT*1 

WRITE (6,700) INLINE 

GO TO 460 

CONTINUE 

IEGFPI=I1EOF+1 | 

GO TO (400,380,370, 380,490, 390} , IEOFP1 
CONTINUE 

CALL HELP (7) 

GO TO 340 

CONTINUE 

WRITE (6,650) 

GO TO 340 

CONTINUE 

IOLIST=NOINPT-ITOTIN 

WRITE (6,600) IOLIST 

IF (IMOOE.NE.1) ICOUNT=ICOUNT+1 
NAMT=MAXINP 

GO TO 410 

CONTINUE 


IF (NOINPT.EQ.ITOTIN) GO TO 410 
ITOTIN=MAXINP-NAMT 

WRITE (6,590) ITOTIN, IROW 

GO TO 350 
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151 
1Se 
153 
154 
155 
156 
15? 
158 
1$S 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
172 
173 
174 
175 
176 
177 
178 
179 
160 
181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
200 


410 


Cc 
902 


420 


430 


“YO 


450 
460 
470 


480 


490 


500 
510 


520 
530 
S40 
550 
560 


570 
580 
590 


600 


610 
620 
630 
640 


CONTINUE 
ILAST=IROW+NAMT-1 


IRERR=0 


00 420 ICOL=IROW, ILAST 


WRITE (6,902) 
FORMAT (°* 
CALL REASON 


IROW, ICOL. ILAST, IMODE, 
ACC’, 4110, 4F10.3) 
(2, ICOL, ]ROW, TEMP (1COL) ,O, IFATAL) 


TEMP (ICOL) 


IRERR=IRERR+I FATAL 


WRITE (6,901) 


TROW, ICOL, IRERR, IFATAL 


CONTINUE 

IF (IRERR.EQ.0) GO TO 44HO 
IF (IMODE.EQ.2) GO TO 430 
WRITE (6,640) IROW 

GO TO 330 
ICOUNT=1COUNT+IRERR 

GO TO 460 

CONTINUE 


DO 450 ICOL=IROW,NMEAS 
N=NOX (TROW, ICOL) 
CORR (N) = TEMP CICOL) 


CONTINUE 

CONTINUE 

CONTINUE 

IF (ICOUNT.EQ.0) ISAMP=1 

IF (IMGOE.NE.1) WRITE (6,680) ICOUNT 

GO TO 490 

WRITE (6,660) 

RETURN 

CONTINUE 

RETURN 

FORMAT (5X, ‘ENTER PAIRS GF MEANS ANO STANOARD OEVIATIGNS--'‘) 

FORMAT (//,° OIAGNOSTICS GF FORMATTED REAO *,/,*° GPER DESC: *,10A4 
if) 

FORMAT {//,1X, ‘CORRELATION/ANTHROPOMETRIC MEASUREMENTS: ‘) 

FORMAT (* MEASUREMENT - ‘,12,° ?---') 

FORMAT (* MEASUREMENT -‘,12) 

FORMAT (5X, "ENTER STANOARD OEVIATION --') 

FORMAT (5X, ‘ENTER ANTHROPOMETRIC CORRELATION MATRIX OATA IN A ROWW 
1ISE UPPER *,/,5X, ‘TRIANGULAR FORM. *,/,5X, "NOTE: ALL DIAGONAL ELEM 
2ENTS MUST EQUAL 1.000") 

FORMAT (SX,4HROW ,12,°C(',12,° ELEMENTS, COLUMNS ‘,1J2,*-',12,°) ‘) 

FORMAT (SX,4HROW ,12,°(',12,° ELEMENT, COLUMN °,12,°') ‘) 

FORMAT (10X,23H»xENTER THE REMAINING ,12,18H ELEMENTS FOR ROW ,12 
1) 

FORMAT (10X, *x»*ACCPTOS1 THE LAST *,1]2,* INPUTS HAVE BEEN IGNORED‘ 
1) 

FORMAT (SX, ‘INPUT MQOE (1=INTERACTIVE; 2=FILE) --*) 

FORMAT (SX, *REAO MOOE (1=FORMATTED;2=UNFORMATTED) -- *) \ 

FORMAT (10X, ‘*»*ACCPTO6O RE-ENTER OATA FOR *,12) | 

FORMAT (10X, ‘xx*ACCPT1300 RE-ENTER DATA FOR ROW ‘, 13) 


oe 





201 
202 
203 
204 
205 
206 
207 
208 
209 


650 
660 
670 


680 
690 
700 
710 


FORMAT 
FORMAT 
FORMAT 


(10X, *xxxACCPT155 INVALIO RESPONSE *) 
(10X, *xx»ACCPT UNEXPECTED EOF ON NIMEAS FOR MEAS °*, 12) 
(10X, ‘x*xREADFLYO30 FILE ANTSTA 1S NOT ATTACHED TO THE PROG 


JRAM*,/,10X, ‘xxxUNFORMATTED READ ON THIS FILE IS ABORTED ---‘//) 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
END 


(10X, *xx»ACCPT3000 TOTAL ERRORS ON INPUT DATA *,15) 
{* MEAN.STD DEV &CORRELATION DATA READ") 

(JOX, ‘x*xACCPT15O0 INVALID DATA AT LINE *,15) 

(10X, *x»x»ACCPTI1200 UNEXPECTED END OF FILE") 


aa 





On Ouwnk why = 


100 


110 


120 


130 


1NO 


150 


160 


170 


180 


190 


200 


210 


220 


SUBROUTINE EDIT 

COMMON/MCSM/CORR (78) , ANTM (2, 12), V¥ (25) , NMEAS, KN, 
SOROOT (78) ,X (25), 1U,NTMEAS, ICOUNT, 
INAME (10) , ITMEAS, NIMEAS, MAXSAM,NTEMP 

COMMON/MISC/IGIN, I00UT, ICARD 

DIMENSION BUFF (3) 

DATA BUFF /3™0./ 

IOIN=5 

CONTINUE 

WRITE (6,350) 

CALL INPT (1,1,NVALSZ,IOPT,XVAL, IEOF) 

JEQFP1=IEOF +1 

GO TO (130,100,110,120,340,120), IEOFP! 

CONTINUE 

CALL HELP (8) 

GO TO 100 

CONTINUE 

WRITE (6,360) 

GO 70 100 

CONTINUE 

IF (JOPT.LT.1.0R.10PT.GT.5) GO TO 120 

GO TO (140,140,240, 330, 340) , 1OPT 

CONTINUE 

ICOL=10PT 

CONTINUE 

MAXINP=e 

NAMT=0 

CONTINUE 

IF (ICOL.EQ.1) WRITE (6,380) 

IF (]COL.EQ.2) WRITE (6,390) 

CONTINUE 

IVALSZ=MAX INP-NAMT 

J=NAMT +1 

CALL INPT (2,1VALSZ,NVALSZ, IVAL. BUFF (J), IEOF) 

JEQFP1=IEOQFP1 

GO TO (220,180,2e00,190,100,210), IEGFP! 

CONTINUE 

WRITE (6,370) 

GO TO 160 

CONTINUE 

WRITE (6,360) 

GO TO 160 

CONTINUE 

IF (ICOL.EQ.1) CALL HELP (11) 

IF (1COL.EQ.2) CALL HELP (12) 

GO 170 160 

CONTINUE 

IOLIST=NVALSZ-IVALSZ 

WRITE (6,400) IOLIST 

CONTINUE 
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230 


240 


250 


260 


270 


260 


290 


300 


310 


320 


IROW=BUFF (1) 

NAMT=NVALSZ+NAMT 

IF (NAMT.GE.MAXINP) GO TO 230 

IF (ICOL.EQ.1) WRITE (6,420) IROW 

IF (ICOL.EQ.2) WRITE (6,430) IRON 

GO TO 170 

CONTINUE 

CALL REASON (1,1COL, 1ROW, BUFF (2) ,0, IFATAL) 
IF (CIFATAL.EQ.1) GO TO 150 

ANTM CICOL, IROW) =BUFF (2) 

GO TO 100 

CONTINUE 

MAXINP23 

NAMT=0 

CONTINUE 

WRITE (6,410) 

CONTINUE 

IVALSZ=MAXINP+NAMT 

J=NAMT +1 

CALL INPT (2, I1VALSZ,NVALSZ,IVAL, BUFF (J), IEOF) 
JEOFP1=JEOF +1! 

GO TO (310,270,290, 280, 100,300), IEOFP! 
CONTINUE 

WRITE (6,370) 

GO TO 250 

CONTINUE 

WRITE (6,360) 

GO 70 250 

CONTINUE 

CALL HELP (13) 

GO 10 250 

CONTINUE 

IOLIST=NVALSZ-IVALSZ 

WRITE (6,400) IOLIST 

GO 70 320 

CONTINUE 

ICOL=BUFF (1) 

]ROW=BUFF (2) 

NAMT=NAMT*NVALSZ 

IF (NAMT.GE.MAXINP) GO TO 320 

IF (NAMT.EQ.1) WRITE (6,440) ICOL 

IF (NAMT.EQ@.2) WRITE (6,450) ICOL,1ROW 
GO TO 260 

CONTINUE 

CALL REASON (3, 1COL,1ROW, BUFF (3) .0, IFATAL) 
IF (IFATAL.EQ.1) GO TO euO 

CALL REASON (2, 1COL, IRON, BUFF (3) ,0, IFATAL) 
IF (IFATAL.EQ.1) GO TO 240 
N=NDX (1 ROW, ICOL) 

CORR (N)} =BUFF (3) 


OFS. 





101 GO TO 100 


102 330 CONTINUE 

103 CALL PRINT 

104 G6 TO 100 

105 34YO CONTINUE 

106 RETURN 

107 350 FORMAT (SX, “EDIT MODE (1-MEANS; 2-STO; 3-ANTH. CORR; Y-PRINT;S-END) -- °) 
108 360 FORMAT (10X, ‘*xxxEQITOWY INVALIO RESPONSE *) 

109 370 FORMAT (10X, ‘xx»EQ01TO9S NO INFORMATION ENCOUNTERED“) 

110 380 FORMAT (5X, ‘INPUT INDEX AND MEAN VALUE *) 

111 390 FORMAT (SX, “INPUT INDEX AND STANDARD DEVIATION*) 

112 YOO FORMAT (10X, *x~xEQIT OWS THE LAST’, J2.° INPUTS HAVE BEEN IGNORED") 
113 4Y1Q FORMAT (5X, ‘INPUT COLUMN . ROW AND CORRELATION*) 

114 4Y20 FORMAT (SX, ‘INPUT MEAN VALUE FOR INDEX ‘*,13) 

115 430 FORMAT (SX, ‘INPUT STANDARD DEVIATION FOR INDEX *,13) 

116 YUOQ FORMAT (5X, “INPUT ROW ANDO CORRELATION FOR COLUMN °,12) 

11? 4YSO FORMAT (5X, ‘INPUT CORRELATION FOR ROW *,12,° COLUMN *, 12) 

118 END 
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100 
110 


120 
130 
140 
150 


160 
170 
160 
190 


SUBROUTINE NAMEFL (1ERROR) 

COMMON/MCSM/CORR (78) , ANTM (2,12) ,V (25) .NMEAS, KN, 
SORCOT (78) ,X (25), 1U,NTMEAS, ICOUNT, 
INAME (10), 1 TMEAS,NIMEAS, MAXSAM,NTEMP 

COMMON/MISC/IOIN, 1O0UT, ICARD 

DIMENSION IBUFF (10), XVAL (10) 

BATA IBUFF/10™4H / 

BATA IYBLK/4UH / 

IERROR=0 

BO 100 J=1,4 

IBUFF (J) =14YBLK 

CONTINUE 

CONTINUE 

IF (10IN.€0.5) WRITE (6,180) 

CALL INPT (4,10,NVALSZ, IBUFF,XVAL. IEOF) 

IF (1EOF.LT.2.0R.1EOF.EQ.5) GO TO 150 

IF (10IN.NE.5) GO TO 140 

IEGFMI=1E0F-1 

GO TG (120,130,170), IEOFM! 

CONTINUE 

CALL HELP (14) 

GO TO 110 

CONTINUE 

WRITE (6, 190) 

G6 TG 110 

CONTINUE 

WRITE (6,190) 

TERROR=1 

CONTINUE 

BO 160 J=1,10 

INAME (J) =IBUFF (J) 

CONTINUE 

RETURN 


FORMAT (5X, “GPERATGR SAMPLE DESCRIPTION (YQ CHAR MAX) ----*) 


FORMAT (10X, *x*xNAMEFLO27 INVALID RESPONSE ‘) 
END 


OF 





SUBROUTINE OPSAMP (ISAMP) 
COMMON/MCSM/CORR (78) , ANTM (2, 12) Vv (25) ,NMEAS,KN, 
j SQROOT (78) ,xX (25) , IU, NTMEAS, ICOUNT, 
2 INAME (10) , ITMEAS,NIMEAS, MAXSAM, NTEMP 


oon ouwuckuwhvo = © 


100 


110 


120 


130 


140 


150 


160 


170 


180 


190 


200 


COMMON/MISC/IGIN, 100UT, ICARD 
DIMENSION IRAN (9) 
DATA IRAN (3), 1RAN (2), TRAN (3), IRAN (4), TRAN (S), IRAN(6) ,ITRAN(7) , 


JIRAN (8) , IRAN (9) /33333333,55555555, 77777777, 222222221, 


UYUNNNNU3, 666666669, 888888889, 999999999, 123456789/ 
IOINSS 
DO 100 K=1,KN 
SQROOT (K) =CORR {K) 
CONTINUE 
CALL SQRMX {(NERROR) 
IF (NERROR.EQ.1) RETURN 
CONTINUE 
WRITE (6, 320) 
CALL INPT (1,1,NVALSZ.NSAMP, XVAL, IEOF) 
IEQFP1=1EOF+1 
GO TO (150,120, 130,140,290, 140), JEOFP! 
CONTINUE 
WRITE (6,330) 
GO TO 110 
CONT INUE 
CALL HELP (15) 
GO TO 110 
CONTINUE 
WRITE (6, 340) 
GO TO 110 
CONTINUE 
IF (NSAMP.GT.O.AND.NSAMP.LE.MAXSAM) GO TO 160 
WRITE (6,350) MAXSAM 
GO TO 110 
CONTINUE 
WRITE (6,300) 
CALL INPT (1,1,NVALSZ,MTYPE, XVAL, IEOF) 
IEOFP121EQF+ 
GO TO (190,160,170,180, 290, 180), IEOFP! 
CONTINUE 
CALL HELP (20) 
Go To 160 
CONTINUE 
WRITE (6, 340) 
GO TO 160 
CONTINUE , 
IF (MTYPE.NE.32.AND.MTYPE.NE.60) GO TO 180 
CONTINUE 
WRITE (6, 310) 
CALL INPT (1,1,NVALSZ,15, XVAL, IEOF) 
IEQFP1=1EQF+1 
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210 


220 


230 


250 


260 


270 


280 


370 


290 


300 
310 
320 
330 
340 
350 
360 


Soa clon etOneceO,230,290,¢30) , IEGFP1 

CONTINUE 

WRITE (6,330) 

GO 70 200 

CONTINUE 

CAE HELP fel) 

GO 10 200 

CONTINUE 

WRITE (6, 340) 

GO T0 200 

CONTINUE 

IF (1S.E0.0) GO TO 250 

trae. 1. 0m.1o.61.9) GE TO 230 

]U=ITRAN (1S) 

GO 16 260 

CONTINUE 

]1U=67111133 

CONTINUE 

REWIND NTEMP 

WRITE (NTEMP) INAME,NMEAS 

DO 280 I=1,NSAMP 

CALL SIMSUB (MTYPE) 

00 270 K=1,NMEAS 

V (K) =V (K) ANTM (2,K) *ANTM(1,K) 

CONTINUE 

WRITE (NTEMP) (V (10) ,10=1,NMEAS) 

CONTINUE 

ISAMP=1 

ENO FILE NTEMP 

REWINO NTEMP 

REWIND NTMEAS 

WRITE (6,370) 

FORMAT ( 10X, ‘x SAMPLE OATA SAVEO ON FILE 9°) 

WRITE (6,360) 

CONTINUE 

RETURN 

FORMAT (5X, ‘ENTER MACHINE WORO SIZE (32 OR 60 BIT) --*) 
FORMAT (5X, ‘ENTER RANDOM NO SEED (0 IF NONE OR 1-9)--") 
FORMAT (5X, ‘INPUT SAMPLE SIZE --‘) 

FORMAT (10X, *x»xxG@PSAMP215 NG INFORMATION ENCOUNTERED *) 
FORMAT (10X, ‘xxxOPSAMP4Y11 INVALIO RESPONSE *‘) 

FORMAT (10X, *x*x»xQPSAMPOOS SAMPLE SIZE MUST BE BETWEEN 1 AND *,13) 
FORMAT (10X, *xx*SAMPLE DATA SAVED ON FILE SAMPLE OR TAPE 3.) 
ENO 


1, 
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100 


110 


120 


130 


140 


150 


160 


170 


180 


SUBROUTINE PRINT 
COMMON/MCSM/CORR (78) , ANTM (2,12) ,V (25) .NMEAS,KN, 
SQROOT (78) ,X (25) , IU,NTMEAS, ICOUNT, 
INAME (10) , ITMEAS, NIMEAS, MAXSAM, NTEMP 
COMMON/MISC/IOIN, 100UT, ICARO 
INTEGER FMT,FMT1 
DIMENSION FMT (6) ,NUM(12) , TEMP (12), INDEX (2) , FMT] (5) 
DATA FMT! /°(1X,"°, ° Oo,  icXkgame Ae,e « XI) */ 
DATA FMT /‘(1X,", ° sO) 5 oe Osa g oy / 
DATA NUM/1H1,1H2,1H3,1H4,.1H5,1H6, 1H7,1H8,1H9,2H10,2H11,2eH12/ 
FMT 1 (2) =NUM (NMEAS) 
ISTORE=0 
]OFLAG=0 
CONTINUE 
WRITE (6,630) 
CALL INPT (1,1.NVALSZ,]OPT, XVAL, IEOF) 
JEGFP1=J]EGF+1 
GO TO (120,110,110.110, 390.110) , JEGFP1 
CONTINUE 
WRITE (6,650) 
GO TO 100 
CONTINUE 
IF (JOPT.LT.1.0R.10PT.GT.4)I GO TO 110 
IF (IGPT.EQ.1) GO TO 130 
IF (IGPT.EQ.2) GO TO 150 
IF (IGPT.EQ.4) GG TO 390 
GO TO 190 
CONTINUE 
WRITE (6,580) 
DO 140 ]ROW=1,NMEAS 
WRITE (6,640) IRGW, (ANTM (JG, IROW) ,JG=1.2) 
CONTINUE 
WRITE (6,620) 
GG 10 100 
CONTINUE 
WRITE (6,590) 
WRITE (6,FMT1) (NUM(I1G),10=1,NMEAS) 
WRITE (6,600) 
DO 180 IRGW=1,NMEAS 
DO 160 ICOL=1,NMEAS 
N=NDX (] ROW, ICOL) 
TEMP (ICGL) =CORR (N) 
CONT INUE 
IF (IROW.GT.1) GO TO 170 
WRITE (6,610) (TEMP (ICOL) , ICOL=] ROW, NMEAS) 
GG TG 180 
FMT (2) =NUM (IRGW-1) 
FMT (4) =NUM (NMEAS*1-1 ROW) 
WRITE (6,FMT) (TEMP CICOL) , ICOL=]ROW, NMEAS) 
CONTINUE 


nO) 





51 
92 
53 
SY 
55 
56 
37 
38 


60 
61 
62 
63 
6u 
65 
66 
67? 
68 
69 
70 
71 
ve 
3 
74 
75 
76 
77 
78 


80 
81 
82 
83 
BY 
85 
86 
87 
88 
89 
90 
91 
92 
aS 
Ou 


96 
97 
98 
Jo 
100 


190 


200 


210 


220 


230 


240 


250 


260 


270 


280 


290 


300 


WRITE (6,620) 

GO TO 100 

CONTINUE 

WRITE (6,430) 

IF (]OFLAG.LT.1) WRITE (6,440) 
]OFLAG=1 

CALL INPT (1,2,NVALSZ, INDEX. XVAL, IEOF) 
JEOFP1=I]EOF+) 

GO TO (230,270,.200.210.100.220), IEOFP! 
CONTINUE 

CALL HELP (19) 

]OFLAG=0 

GO TO 190 

CONTINUE 

WRITE (6,420) 

GO TO 190 

CONTINUE 

IOLIST=NVALSZ-2 

WRITE (6,450) IOLIST 

CONTINUE 

IF (INOEX (1).GT.0) GO TO 240 

WRITE (6,480) 

GO TO 190 

CONTINUE 

IF (CISTORE.GT.O.ANO.INOEX (1) .GT.ISTORE) GO TO 250 
NRSKIP=INOEX (1) -1 

GO TO 260 

CONTINUE 

NRSKIP=INDEX (1) -ISTORE-1 

CONTINUE 

ISTART=INOEX (1) 

IF (NVALSZ.EQ.1) GO TO 280 

IF (NVALSZ.GE.2) GO TO 290 

CONTINUE 

ISTORE=0 

ISTART=1 

IENO=999 

NRSKIP=0 

GO TO 320 

CONTINUE 

TENO=ISTART 

GO T0 310 

CONTINUE 

IF (INOEX (1) .LE. INDEX (2)} GO TO 300 
WRITE (6,490) 

GO TO 190 

CONTINUE 

IENO=INDEX (2) 

CONTINUE 

IF (INOEX (1) .GT.ISTORE.ANO.ISTORE.NE.O) GO TO 330 


TOL 





101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
lle 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
131 
ise 
133 
134 
135 
136 
137 
138 
139 
140 
141 
ie 
143 
1u4 
145 
146 
147 
148 
149 
150 


320 


330 


34D 
350 


36D 


370 


380 


390 


400 


410 


420 
43D 
“UD 


CONTINUE 
REWIND NTEMP 


READ 


CONTINUE 


IF 


(NRSKIP.LT.1) 


(NTEMP,ENO=410) INAME,NMEAS 


GO T0 350 


DO 34D I=1,NRSKIP 

READ (NTEMP,END=40D) (V{K) ,K=1,NMEAS) 
CONTINUE 
CONTINUE 
WRITE (6,530) (INAME (10), 10=1, 10D) 

(ISTART.EQ.1ENO) WRITE (6,540) IENO 
(ISTART.NE.TEND.AND.ITEND.NE.999) WRITE (6,550) ISTART, IENO 
WRITE (6,500) 
WRITE (6,FMT1) (NUM (10), 10=1,NMEAS) 
WRITE (6,600) 
ISKIP=0 
NFOUNO=0 

DO 360 I=ISTART,IEND 


IF 
IF 


IF 


(ISKIP.EQ. 1) 


ISKIP=D0 
(ISTART.NE.TENO.AND.I.NE.ITENO.ANO. (1/5) *5.£0.1) ISKIP=1 


IF 


REAO 


NF OUNO=NFOUNO*! 
ISTORE=1 

WRITE (6,510) (V(K) ,K=1,NMEAS) 
CONTINUE 
GO TO 380 
CONTINUE 


IF 


(TEND.NE.999. 


ISTORE=1 


IF 


(NFOUND.LT.1) 


CONTINUE 
WRITE (6,520) 
GO T0 190 
CONTINUE 
REWINO NTEMP 
RE TURN 
CONTINUE 
IOLIST=I-2*I1STORE 
WRITE (6,560) 
GO 170 190 
CONTINUE 
WRITE (6,570) 
GO 10 100 


FORMAT 
FORMAT 
FORMAT 


1‘FOR INDIVIDUAL 
2A SEQUENCE OF OPERATORS (ENTER RANGE OF OPERATCGR NUMBERS) *,/7,10X, ° 


3ENO 


(ENTER $) *) 


(NTEMP,ENO= 


WRITE (6,46D) 


37D) (V(K),K=1,NMEARS) 


AND. ISTART.NE.IEND) WRITE (6,470) IEND 


WRITE (6,560) 


(10X, *xxxPRINTI3S5 INVALIO RESPONSE ‘) 
(SX, ‘OUTPUT MODE: --*) 
(10X, ‘FOR ALL OPERATOR DATA (ENTER CARRIAGE RETURN} *, 7, 10X, 


OPERATOR DATA (ENTER GPERATOR NUMBER) ‘,/7,10X, ‘FOR 


OZ 





io! 
152 
153 
154 
155 
156 
157 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
173 
172 
173 
174 
175 
176 
17? 
178 
179 


4¥S0 
460 
470 
480 


490 


500 
510 
520 
530 
S40 
550 
560 


570 
580 
590 
600 
610 
620 
630 


640 
650 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
10 ZERC 
FORMAT 
1X*) 

FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 


(10X, *x«*xPRINT119 THE LAST‘, 14, “INPUTS HAVE BEEN IGNORED*) 
(1H ) 

(10X, *xxxPRINT211 ERROR I§ UPPER INDEX --*, 14) 

(10X, *«x*PRINTI20 INDEX NO. LESS THAN GR EQUAL’, /, ‘xT 

1S NOT VALIO*) 

(10X, *xx»PRINTI60 FIRST INDEX MUST BE LESS THAN SECOND INDE 


(27X, ‘OPERATGR SAMPLE*,/, 71 (1H-)) 

(11F6.2,F5.e2) 

(71 (1H-)) 

(1X, 10A4, 10X,AY) 

(10X, ‘OPERATGR NO*, 14) 

(10X, “OPERATGR NOS*,14,* --°, 14) 

(10X, "xx*PRINT1S00 INDEX NUMBER EXCEEDS RECORDS IN OPERATOR 


1 FILE *) 


FORMAT 


(10X, *xx*PRINTI600 FILE SAMPLE 15 NOT ATTACHED TO THE PROGR 


1AM*,/7,10X, *xxxGPERATQR GUTPUT REQUEST ABORTED ---"*,//) 


FORMAT 


(7,21X, ‘MEANS AND STANDARD OEVIATIONS*,/,1X,70(1H-),/,2X, “M 


1EASUREMENT*,6X, “MEANS *,6X, “STO OEV*,/,1X, 71 (1H-) ) 


FORMAT 


(/,7X, ‘CORRELATION MATRIX FOR ANTHROPOMETRIC MEASUREMENTS", 


1/,1X, 70 (1H-)) 


FORMAT 
FORMAT 
FORMAT 
FORMAT 
10)--"*) 
FORMAT 
FORMAT 
ENO 


(1X, 70 (1H-)) 

(1X, 12F 6.3) 

(1X, 70 (1H-) ) 

(SX, "PRINT MODE (1-MEAN/STO;2-CORR MATRX; *. °3-QPER MEAS: 4-EN 


(7X, 12,5X,2 (SX, F7.2)) 
(10X, *x«xPRINT103 INVALIO RESPONSE *) 


ines 
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100 


110 


SUBROUTINE RANDe (RNORM,MTYPE, QUAN, IJ) 

COMMON/MCSM/CORA (78) ANTM (2, 12) ,V (25) ,. NMEAS, KN, 
SQROOT (78) .X (25) , IU,NTMEAS, ICOUNT, 
INAME (10) , ITMEAS,NIMEAS, MAXSAM, NTEMP 

DIMENSIGN QUAN (2) 

COMMON/MISC/IOIN, 1GGUT, ICARD 

IJ=]J+1 

IF (1J.EQ.2) GO TO 110 

CONTINUE 

X1=2.0xUNFORM (MTYPE) -1.0 

Y1=2.0™UNFGRM (MTYPE) -1.0 

S=X1eX1+VixV1 

IF (S.GE.1.0) GO TO 100 

S=SQRAT (-2. OxALOG (S$) /S) 

QUAN (1) =X1™5 

QUAN (2) =Y1%5 

CONTINUE 

RNORM=QUAN (1 J) 

IF (I1J.EQ.2) IJ=s0 

RETURN 

END 


104 
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901 


100 


110 


120 


130 


140 


150 


160 


SUBROUTINE REASON (ICALL, ICOL, IROW, TSTOAR!, TSTORe, IFATAL) 
COMMON/MCSM/CORR (78) , ANTM (2, 12) ,V (25) ,NMEAS,KN, 
1 SQROOT (78) ,X (25) ,1U,NTMEAS, ICOUNT, 


2 INAME (10), ITMEAS, NIMEAS, MAXSAM, NTEMP 


COMMON/MISC/IGIN, 100UT, ICARO 

IFATAL=0 

IAFLAG=0 

WRITE (6,901) ICALL,ICOL, IROW, IFATAL, TSTOA], TSTORS 
FORMAT (* REASON *, 4110, YFIO. 3) 

IF (ICALL.LT.O) IAFLAG=1 

ICALL=IABS (ICALL) 

IF (ICALL.GT.1) GO TO 110 

IWRITE=] 

IF (CIROW.LT.1.0R.IROW.GT.NMEAS) GO TO 120 
IF (1COL.EQ.2) GO TA 100 

RETURN 

CONT INUE 

IWRITEs2 

IF (TSTGR2.LT.O.) GO TO 120 

IWRITE=3 

IF (TSTORI.LE.TSTOR2) GO TO 120 

RETURN 

CONTINUE 

IWRITE=1 

IF (]ROW.LE.O.0R.1ROW.GT.NMEAS) GO TO 130 
IWRITE=2 

IF (I1COL.LE.0.QOR.1COL.GT.NMEAS) GO TG 130 
IWRITE=3 

IF (IROW.GT.ICOL) GO TO 130 

IF (I1COL.EQ.3) RETURN 

IWRITE=4 

IF (ICOL.EQ. IROW.AND. TSTORI.NE.1.000) GO TA 130 
IWRITE=5 

IF (TSTORI.LT.-1.000) GO TO 130 

IWRITE<6 

IF (TSTOR1.GT.1.000) GO TO 130 

RETURN 

CONTINUE 

IF CIWRITE.NE.1.AND.JQIN.NE.5) WRITE (6,280) IROW 
GO TO (190,140,150), IWRITE 

CONTINUE 

G6 TAO (190,200,210,160,170,180), IWRITE 
CONTINUE 

WRITE (6,230) 

GO 10 220 

CONTINUE 

WRITE (6,240) 

GO 10 220 

CONTINUE 

WRITE (6,250) IRGW, ICOL 
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50 
51 


53 
54 


56 
57 
58 


60 
6 
62 
63 
64 


66 
67 
68 
69 
70 
71 
ie 
73 
7u 
75 


77? 
78 
79 
60 
81 
82 
83 
Bu 


170 


180 


190 


200 


210 


220 


230 
240 
250 
260 
270 


280 
290 


300 


310 


GO TG 220 

CONTINUE 

WRITE (6,260) IROW, ICOL 

GO T@ 220 

CONTINUE 

WRITE (6,270) IROW, ICOL 

GO TO 220 

CONTINUE 

WRITE (6,290) IROE,NMEAS 
GO TO 220 

CONTINUE 

WRITE (6,300) ICOL,NMEAS 
GO TO 220 

CONTINUE 

WRITE (6,310) IROW, ICOL 

GO TO 220 

CONTINUE 

IFATALS1 

RETURN 

FORMAT (10X, °xxxREASONI 150 
FORMAT {(10X, °xx*REASON1I200 
FORMAT {(10X, °xxxREASONI2S0O 
I1COLUMN °,12,° 1S NOT EQUAL 
FORMAT (10X, *xx»xREASON1300 
iumxeIS LESS THAN —1.0°) 
FORMAT (10X, °x**REASON1350 
fux]S GREATER THAN +1.0°) 
FORMAT {(8X, 12) 

FORMAT (10X, °***REASON1SO0 
1 THAN °,12) 

FORMAT (10X, °~x*xREASONI600 
ITER THAN °,1J2) 


STANOARD OEVIATIONS LESS THAN ZERO’) 
STO. GREATER THAN MEAN VALUE. °) 
DIAGONAL ELEMENT IN ROW *,12,/7,10X, ’x 
105)20)7) 
ENTRY IN COLUMN °,J]2,° ROW *,12,/,10X, ‘ 


ENTRY IN COLUMN °,12,° ROW *,12,/7,10X,° 


ROW *,12,° IS LESS THAN ZERG OR GREATER 


COLUMN 


velo. 9 


1S LESS THAN ZERO OR GREA 


FORMAT (10X, ‘xx*REASONI 700 ROW °,12, ‘COLUMN °,]2,° 1S GUTSIDE OF 
1UPPER*,/,10X, ‘xx*xTRIANGULAR REGION OF CORRELATION MATRIX. °) 


ENO 


106 





won on u£& wh = 


100 


SUBROUTINE SAVE 
COMMON/MCSM/CORR (78) , ANTM (2, 12) .V (25) ,NMEAS,KN, 
i SQROOT (78) .X (25) , 1U,NTMEAS, ICOUNT, 
2 INAME (10), ITHEAS,NIMEAS,MAXSAM, NTEMP 
COMMON/MISC/IOIN, IOOUT, ICARO 

REWIND NTMEAS 

WRITE (NTMEAS) NMEAS, KN, INAME 

WRITE (NTMEAS) ( (ANTM (10, JO), 10=1,2) ,JO0=1,NMEAS) 
WRITE (NTMEAS) (CORR (10), 10=1,KN) 

REWIND NTMEAS 

WRITE (6,100) 

RETURN 

FORMAT (* MEASEUREMENT DATA SAVED ON FILE ANTSTA GR TAPE 2°) 
END 


Od, 
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100 


110 


SUBROUTINE SIMSUB (MTYPE) 
COMMON/MCSM/CORR (78) , ANTM (2, 12), V (25) , NMEAS, KN, 
1 SQROOT (78) ,X (25) , 1U,NTMEAS, ICOUNT, 


2 INAME (10), ITMEAS, NIMEAS, MAXSAM,NTEMP 


COMMON/MISC/IGIN, 100UT, ICARD 
DIMENSION Y (45) 

IJ=0 

00 100 1I=1,NMEAS 

V (I) =0.0 

CALL RANDe (RNORM,MTYPE. QUAN, IJ) 
Y (1) =RNORM 

00 110 JI=1,NMEAS 

00 110 J=1,1 

K=NDX (J, 1) 

V (1) =V 01) +¥ (3) *SQROOT (K) 
RETURN 

END 
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100 


110 


120 
130 


140 


150 


160 


170 


SUBROUTINE SQRMX (NERROR) 

COMMON/MCSM/CORR (78) , ANTM (2,12) .V (25) ,NMEAS, KN, 
SQROOT (78) ,.X (25) , IU,NTMEARS, ICOUNT, 

INAME (10), 1 TMEAS,NIMEAS, MAXSAM, NTEMP 

COMMON/MISC/IOIN, IJGOUT, ICARD 

DIMENSION 2X {7) 

INTEGER 2X 

NERROR=0 

DO 120 I=1,NMEAS 

ZX (4) =] 

N=NDX (1,17) 

IF (SQROOT(N) .LT.O) GO TO 130 

SO=SQRT (SQROOGT (N) ) 

DO 100 J=I,NMEAS 

N=NDX (J, J) 

SQROOT (N) =SQRCGGT (N) /SQ 

CONTINUE 

IF (].£Q.NMEAS) NERROR=0 

IF (1.EQ.NMEAS) RETURN 

li=]+1 

ZX (5) 3] 1 

DO 110 Ie=11,NMEAS 

DO 110 Je=12,NMEARS 

N=NOX (Ie, Ja) 

L=NDX (1,12) 

K=NDX (I, Je) 

SQROOT (N) =SQROOT (N) -SQROOT (L) *SQROOT (K) 

ZX (6) =12 

ZX (7) =Je 

CONTINUE 

ZX (1) =N 

ZX (2) st 

ZX (3) =K 

CONTINUE 

CONTINUE 

NERROR=! 

WRITE (6,180) 

CONTINUE 

WRITE (6,190) 

CALL INPT (4,1,NVALSZ, IANS, XVAL, IEOF) 

TEQFP1I=1EOF+1 

GO TO (170,160,150,170,160,170), IEOFP! 

CONTINUE 

WRITE (6,200) 

WRITE (6,210) (ZX (1),1=1,7),SQROGT (N) 

CONTINUE 

RETURN 

CONTINUE 

WRITE (6,220) 

GO TO 140 


ee 





180 FORMAT (10X, ‘*xx*SQRMX051 ANTHROPOMETRIC CORRELATION VALUES ARE INV 
1AL10O*) 

190 FORMAT (10X, ‘-- ENTER CARRIAGE RETURN TO CONTINUE PROCESSING*,/, 10 
1X,*-- ENTER ? TO OISPLAY SORMX VARIABLES. *) 

200 FORMAT (10X, ‘xx*SQRMXO62 AN INDEFINITE QUANTITY HAS BEEN OETECTED 
1BY*,/,10X, *x»~SQRMX. VARIABLE LIST WITHIN SQRMX:*,/7,10X* N E 
2 K ] 1 le J2 SQROOT*,/,10X,44 (1H-)) 

210 FORMAT (8X,4(3X,12) ,.4X,3 (12,3X),F6.3,/,10X, 44 (1H-) ,//) 

220 FORMAT (10X, *xx*xSQRMXO60 INVALID RESPONSE x» *) 
END 


LAO 
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100 


FUNCTIGN NDX (1, J) 

IF (1.GT.J) GG TG 100 
NOX=1] + (JuJ-J) /2 
RETURN 

NDX=J+ (IuT~-1) /2 
RETURN 

END 
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110 


FUNCTION UNFORM (MTYPE) 

COMMON /MCSM/ CORR (78), ANTM (2,12) ,V (25) ,NMEAS, KN, 
1SQROOT (78), X (25) , ]U,NTMEAS, ICOUNT, 
2INAME (10), 1 TMEAS, NIMEAS, MAXSAM,NTEMP 

COMMON/MISC/ICIN, 1OOUT, ICARO 

DATA 1X/16777213/,C/281474976710655./ 

IF (MTYPE.EQ.32) GO TO 100 

1U=]Ux1X 

UNFORM=I1U/C 

RETURN 

1U=1Ux65539 

IF (IU.GE.0) GO TO 110 

1U=1U+21U7NB36N7+1 

UNF ORM=1]U 

UNF ORM=UNF ORM™. 465661 3E-09 

RETURN 

END 


ee 
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SUBROUTINE HELP 
RETURN 
END 


(1) 


113 


rf) 





$9939 ] 


100000 BLOCK DATA 

100001 COMMON/MCSM/CORR (78) , ANTM (2,12) ,V (25) ,NMEAS,KN, 
100002 1 SQROOT (78) ,X (25) ,1U,NTMEAS, ICOUNT, 

100003 2 INAME (10), ITMEAS,NIMEAS, MAXSAM, NTEMP 

100004 COMMON/MISC/IOIN, 100UT, ICARD 

100005 DATA I10IN, 1OOUT,NTEMP,NTMEAS, NIMEAS/S,6,3,4, 7/ 
100006 DATA MAXSAM/400/ 

100007 DATA NTMEAS,NIMEAS/2,4/ 

100008 1, ICOUNT/O/ 

100009 DATA ANTM/24%0./,CORR/78™0./, INAME/UxYH / 
100010 END 

100011 1 
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SUBROUTINE PUNRT | 
COMMON/MCSM/CORR (78) , QNTM (2, 12), ¥ (25) ,NMEAS, KN, 
1 SQROOT (78) ,X (25) , 1U,NTMEAS, ICOUNT, 
2 INAME (10), 1TMEAS, NIMEAS, MAXSAM, NTEMP 


oor ouk wr — © 


27) 


370 


27e 


140 


141 


COMMON/MISC/IOIN, 100UT, ICARD 
REWIND NTEMP 
REWIND 9 


CONTINUE 
READ (NTEMP,END=140) INAME,NMEAS 
WRITE (9,271) INAME,NMEAS 


CONTINUE 

READ (NTEMP,END=370) (YK) ,K=1,NMEAS) 
WRITE (9,271) (V(K) ,K=1,NMEAS) 

FORMAT (20AY) 

GO TO 4Y 

CONTINUE 

ENDFILE 9 

WRITE (6,272) 

FORMAT (10X, ‘xx™~ DATA SAVED ON FILE 9°) 
REWIND NTEMP 

RETURN 

CONTINUE 

WRITE (6,141) 


FORMAT (10X, °* ERROR DATA SET DOES NOT EXIST") 


RETURN 
END 
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